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Pies ReS C1 
The M&ssbauer effect has been used to investipate electronic 
shielding by closed electron shells in salts of trivalent thulium, by 
Se@easuring the temperature dependence of the nuclear quadrupole 


ans, 


Splitting of the 8.42 keV pamma transition in Tm The nuclear 
Guadrupole interaction was studied for igor ions in thulium ethyl 
sulfate, thulium oxide and thulium trifluoride within a temperature 
meee from 9. 6° K to 1970° K. The interpretation of the experimeqses 
mmein terms of the contributions of distorted closed electron shells 
to the quadrupole interaction yields values for electronic shielding 
meetors. The results lead to amounts of 10% or Jess for the atomic 
Seermpneimer factor Roy: The experiments also reveal substantial 
mmmercding of the 4f electrons from the crystal electric field, ex- 
pressed by the shielding factor v5. Values of 250 and 128 are ob- 


famed for the ratio (l-y_)/(l-—c,) for thulium ethyl sulfate and thu- 
ee y 


>) 


lium oxide respectively, where y is the lattice Sternheimer factor. 
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a 
Pee i wo puUCc TION 

The technique of recoilless nuclear resonance absorption of 
gamma radiation, the so-called Missbauer effect (1), has been em- 
ployed in numerous experiments in recent years (2), (3). By binding 
a radioactive nucleus ina crystal Jattice the emitted gamma radiation 
will, under certain conditions (2), have essentially the natural line 
width as determined by the Heisenberg uncertainty relation and an 
energy exactly equal to the excitation energy of the nucleus. Ifa 
nucleus of the same isotope which is in its ground state is also bound 
ina lattice, there is a large probability for nuclear resonance ab- 
sorption of the gamma radiation. A distinct advantage of this tech- 
nique is the inherently high energy resolution that is available. For 
example in the experiments to be described here, the resolution is 
One part in oer Energy resolutions of this order make it possible 
to study nuclear properties as well as solid state effects in the crys- 
tals that are used to bind the nuclei. We make use of the MUssbauer 
effect here to study the nuclear hyperfine interactions in salts of 
rare earths, specifically thulium salts. 

Measurements of the nuclear quadrupole interaction in salts of the 
rare earth elements yield information on the quadrupole moments of 


the relevant nuclear states and onthe electric field gradients which 


(1) R.L. Mussbauer, Z. Physik 151, 124 (1958); Naturwissen- 
pemcnten 45,. 590 (1958), 4. Naturforsch. Idaye2i (1959) 





(2) See for instance H. Frauenfelder, The MUussbauer Effect, 


(W.A. Benjamin Inc., New York, 1962) 


(3) The Proceedings of the Third International Conference on the 
Mussbauer Effect appear in Rev. Mod. Phys. 36, 333-504 (1964) 





ioe 
exist in the salts at the nuclear sites. The extraction of the compo- 
nents of the electric field gradient tensor from such measurements 
is rather straightforward if the values of the nuclear quadrupole mo- 
ments have been obtained by other methods such as Coulomb excita- 
tion techniques. Onthe other hand the determination of nuclear mo- 
ments of rare earth nuclei by measurements of the nuclear quadru- 
pole interaction is rather involved since this requires a calculation 
m@etne components of the electric field pradient tensor at the nuclear 
sites. A calculation of the electric field gradients for salts of the 
rare earths can be performed at present only with limited accuracy. 
Uncertainties in excess of 30% are typical. It therefore appears 
that measurements of the nuclear quadrupole interaction in solids of 
the rare earths are at present of more importance for studies of the 
sources of the electric field gradients than for determination of 
nuclear quadrupole moments. 

The electric field gradient at the nuclear site of a certain ion 
Originates from a number of different sources. Major sources are 
distortions of the electronic shells of the ion. These distortions 
result from the interactions of the electrons of the ion with the crystal 
electric field (CEF) produced by the surrounding ions in the lattice, 
provided the arrangement of the surrounding ions reflects a point 
Synmimetry lower than cubic. The field gradient at the nuclear site 
results not only from the distorted partially filled 4f electron shell 
of the rare earth ion, but also from distorted closed electron shells. 
These distortions of the closed electron shells of the ion constitute a 


major source of uncertainty in caleulations of the electric field 





oe 

Peacdient at the nuclear site. The deviations of the closed shells from 
spherical symmetry (electric multipole polarization) usually lead to 
substantial reduction or enhancement (shielding or antishielding) of 
the electric field gradient at the nuclear site. Sternheimer (4), (5) 
was first to emphasize the importance of magnetic dipole and electric 
quadrupole polarizations of closed shells and pioneered in calculating 
the contributions of closed shell deformations to the nuclear hyper- 
fine interactions. 

The nuclear quadrupole interaction depends strongly on the elec- 
Meonic State of the ion. The electronic states which arise when a 
rare earth ion is incorporated ina crystal lattice are basically 
@aused by the interaction of the CEF and the electrons in the partially 
meeed 4f electron shell, but the splittings of these electronic levels 
are also strongly influenced by distortions of the closed electron 
shells (6),(7),(8). In order to account for the modification of the 
CEF splitting which results from electronic shielding, one has to 
consider the quadrupole moment as well as higher multipole mo- 


ments induced in the closed shells. 


(4) Pei Sternneimenr, Phys. Rev. 80;-102 (1950), [05m woe 
moo); ROM, Sternheimer and H.M. Foley, ibid; T02Z aor 
(1956); H.M. Foley, R.M. Sternheimer and D. Tyko, ibid. 


93, 734 (1954) 





(5) Real, Sternhneimer, Phys. Rev. 84,244 (19> Ten 750, 736 (1954) 
(6) Det eedimonds, Phys. Rev. Letters 10; 129 (1963) 


(7) R.G. Barnes, E. Kankeleit, R. L. MUssbauer and J.M. 
Poindexter, Phys. Rev. Letters 11, 253 (1963) 


(8) J. Blok and D.A. Shirley (private communication) 
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Rare earth ions exhibit CEF splittings which are usually very 
much smaller than similar splittings observed for ions of the iron 
transitions elements. Inthe irontransition series, the partially 
meeed 3d electron shell is fully exposed to the CEF produced by 
Surrounding ions, resulting in large CEF level splittings. The rela- 
tively small CEF level splittings observed for rare earth ions, 
which typically are of the order of a few hundred cm" ', probably 
arise because of large shielding effects resulting from the Sct 
electronic shells which surround the partially filled 4f shell. 

feeesent theoretical predictions of the influence of electronic 
shielding upon the CEF level splitting of rare earth electronic levels 
diverge. Burns (9) concluded that electronic shielding in the rare 
earth ions is of little importance and that the difference between the 
CEF level splittings in the iron scries and those in the rare earth 
series cannot be attributed to electronic shielding of the 4f electrons 
Meem the CEF by outer closed electron shells. In contrast, Lenander 
and Wong (10), Ray (11) and Watson and Freeman (12) conclude that 
electronic shielding plays a significant role in rare earth CEF level 
splittings. 

Quantitative estimates of actual shielding effects are hampered 


by the lack of sufficiently accurate atomic wave functions for rare 


(9) ee burns, Ohys, Rev. 128, 2121 (1962) 
me C.J. Lenander and E. Y. Wong, J. Chem. Phys. 38, 2750 (1963) 
Smee. i. Ray, Proc. Phys. Soc. 82, 47 (1963) 


(12) R.E. Watson and A.J. Freeman, Phys. Rev. 133, A1571(1964) 
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earth ions. Inadequate knowledge of the contributions of the core 
electrons is a primary source of uncertainty in our present under- 
standing of hyperfine interactions in rare carth (as well as in other) 
elements. Direct measurements of the influence of electronic 
ee lding upon the nuclear hyperfine interactions and upon the CEF 
splittings of electronic levels therefore are highly desirable. 

ieais paper demonstrates the use of the technique of recoilless 
nuclear resonance absorption of gamma radiation as a means to ob- 
Eeiantormation on electronic shielding effects in rare earth isotopes. 
The procedure introduced here consists of combining measurements 
mmene temperature dependent nuclear quadrupole interaction (per- 
formed by using the technique of recoilless resonance.absorption) 
with measurements of the CEF level splittings (performed by using 
optical techniques). Specifically we report on determinations of the 
relevant electronic shielding factors for trivalent thulium based 
upon our gamma-absorption measurements of the nuclear quadrupole 


69 


imieeraction of Tent in thulium ethyl sulfate (13) and thulium oxide 


and on optical measurements of CEF levels by Wong and Richman 
(14), Gruber and Krupke (15), and Gruber et al. (16). 


Tm |°9 


appeared to be an isotope particularly suited for studies 


(13) A preliminary report of part of this work appeared eclsewhere 


(7). 
me) &.Y. Wong andl. Richman, J.Chem. Phys. 34, 1182 (1961) 
(15) J.B. Gruber and W.F. Krupke, to be published 


(16) J.B. Gruber, W.F. Krupke and J. M. Poindexter, to be 
published 
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of electronic shielding, for the following major reasons: 

(1) The low energy of the 8.4 keV transition used results ina high 
Debye-Waller factor (recoil-free fraction) even at very high tempera- 
tures, thus permitting a measurement of the quadrupole interaction 
within an unusually wide temperature range. 

(2) The separation of the excited levels belonging to the ground 
multiplet of thulium (L = 6; S = 1) is rather large, with the first 
excited level CH,) some 5600 ae above the ground term Hg). 
Thus the existence of the higher levels of the ground multiplet is of 
minor concern for the interpretation of our data in thulium, in con- 
trast to the situation prevailing in the case of some other rare earth 
ions. 

(3) The spin of the nuclear ground state (I = 1/2) and of the 8.4 keV 
excited state (I = 3/2) is rather low, resulting in a small number of 
quadrupole hyperfine components of the gamma lines which are easily 
resolvable. 

(4) The nuclear collective model applies well to Pra ide thus per- 
mitting a rather reliable semi-theoretica] estimate of the nuclear 
quadrupole moment of the 8.4 keV _ state. 


169 


(5) The relative abundance of Tm is 100%. 
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eee ct vol Al ELECTRIC FIELD (CEF) INTERACTIONS 

A rare earth ion interacts ina salt with the CEF produced by all 
the ions which surround its position in the lattice. The dominant 
effect is the interaction of the CEF with the electrons in the partially 
filled 4f-shell. This interaction is weak compared to the spin-orbit 
interaction, in contrast with the situation prevailing in the case of 
iron-transition elements. Asa result, the total angular momentum 
J remains a good quantum number for rare earth ions bound in 
crystals. The effect of the CEF then essentially is a partial or com- 
plete removal of the 2J + 1 fold spatial degeneracy of the orientation 
Smee which exists ina free ion. The actual number of electronic CEF 
levels depends on the symmetry of the field, while the level spacing 
depends on the strength of the interactions between the CEF and the 
4f electrons. The situation is illustrated in Fig. la. 

The potential energy describing the interaction between the CEF 
and a negative charge at position (r, 3, ¢ ) within the ion centered at 
the origin may be represented in good approximation by the following 


expansion, not including shielding from closed shells: 


vem 
rm. Sas 
-eV(r,U,P)= > r. Nile Nye 5 7) (1) 
n n 
n m=-n 
if one assumes that there is no overlap between the charge distri- 
: - m 

Meetons of different ions. In Eq. (1) the A represent lattice sum. 
Meet point charges and effective multipole moments in the surround= 


.; m : 
ing ions. The relevant functions Qo. , which are linear combinations 





ue. la: 


Pie 


i: 


os 

Schematic of the atomic level splitting of a rare earth 
ion inthe CEF. Fora nuclear spin I = 3/2 the nuclear 
quadrupole interaction splits each CEF level into a 
doublet, which is the case illustrated. Typical overall 
CEF splittings are of the order of oe eV, while typical 
6 


quadrupole hyperfine splittings are of the order of 10° 


eV. 


Schematic of the nuclear quadrupole splitting of the 

8.4 keV transition in ie The temperature depenaens 
level splitting ae: which is typically of the order of 
Non. eV, is the average of the hyperfine splittings of 


Fig. la, weighted according to their Boltzman factors. 
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-1]U= 
- a ere 


Oe e0| Lad : . ate 
pherical harmonics Eo and ve , are defined as follows (17): 


Be Cen) (con a) 


Bet tT) m eek) m cos m? 
@ =(2) (m!) -——— PF (cos vw) . m7? 0 
n (efi) n : 
sin mY 
where P_ and ies are Legendre polynomials and associated Legendre 


functions, respectively. In particular, we obtain for n = 2: 
Srcos- a7 = | (2a) 
sin? #cos2¥ (2b) 


05°= sin? @ sin 2 7 (2c) 

Specifically, the Hamiltonian describing the interaction between 
the CEF and the electrons in the partially filled 4f shell of rare earth 

ions, micluding the eifect of shielding via the closed electron sheiis 


of the central ion is given by: 





m n a 
A PR SiO) PLA) 


(4f) = 
CER a *. 


nm 


(17) The normalization of the functions Oo (uv, Y) is arbitrary; 
the choice adopted here is the one most commonly used 





el Oe 


n : ; 
Sie terms proportional to r, describe the potential energy due te 


k 
the direct interaction of the CEF with the k-th electron in the 4f-shell 
while the terms proportional to Sify) describe the additional poten- 
tial energy arising from a deformation of the closed electron shells. 
The interaction described by the Hamiltonian in Eq.(3) splits the 
electronic ground state of the free ion, characterized by total angu- 
lar momentum J, into a number of CEF levels. We shall assume 
in calculating these CEF levels that the angular and radial parts of 
the free ion wave functions can be factorized and that higher terms 
with different J values can be neglected. Under these circumstances 
we are dealing with a manifold of states belonging to the same J and 
fereethnen convenient to replace the angular operators occurring in 


the Hamiltonian, Eq.(3), by equivalent operators (18). The relevant 


matrix elements then are of the form 
- m n fo 
pS GO Cea On (JJ S,iI.m}? (4) 


where Cy ee\to: CTD. (5) 


(18) MW tin otevens, Proc. Phys. Soc. Aods, 209 (1957), 
R.J. Elliott and K.W.H. Stevens, Proc. Roy. Soc. A218, 
Po 190 5e), 0.8 Diliott, B, RR. Judd and- WA. ekvumeiimae 
pord. 240,909 (19577, R. Orbach, ibid. A264, 458 (1961) 








ml? < 





r= Oye | S00 Uae > [<= 71 (©) 
Cas = Wye | al Ugg? a 
and vy = sa py yo “lor. “1h Z2, ae 


In these expressions U, is the radial part of the electronic wave 


f 
functions for the 4f-shell. The functions one oie Jo) are opera- 
tor equivalents; those relevant for this work are listed in Table I. 
The expressions J | v 5 2 are reduced matrix elements (19), 
which for the more general case of intermediate coupling are avail- 
able for Tel in the literature (14)-(17), (20). 

iis in principle possible to calculate the parameters ae and 
oe p.. but difficult in practice. Difficulties are in the evaiuation 
of the "lattice sums" ae because of a lack of sufficient knowledge of 
the ionic position coordinates and their temperature dependence as 
well as of the values of moments in the surrounding ions (21). The 
evaluation of the radial integrals Cie which are the expectation 
values of r for the 4f shell modified by contributions from closed 
shells to the electric multipole fields at the 4f electron positions, is 


hampered by the lack of knowledge of sufficiently accurate atomic 


(19) Panevudde or roc. Roy. soc.) Ag4l, 414 (193) 


(20) Pea wariber ancl. G. Conway, J. Chem. Phys. o2amocr 
(1960) 


(21) Mi trstinutchings and Di: K. Ray, Proc’ Phys. Soc. 8), 605 
(1903) 
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ive functions for bound rare earth ions. For these reasons it 1s 


therefore preferable to introduce the 'CEF parameters'! 


to be determined by experiment. The point symmetry of the central 
ion drastically reduces the number of CEF parameters (22). In the 
case of rare earth ions only the terms with n= 2,4, 6 need to be 
considered, with the effects of n= 1, 3, 5 being neghyible in most 
cases (22). 


The wave functions G of the V-th CEF level will be taken asa 


linear combination of eigenvectors of the total angular momentum J , 
7 Cnet 


(9) 


The expansion coefficients ay and the energy eigenvalues Ey, 


follow from the diagonalization of the interaction matrix ae 
J, aa 


Eq. (4). 





Nac) A compilation of the relevant values n and m for various 
crystal symmetries was given by J. L. Prather, NBS 
Monograph 19 (1961) 


7 








Be 
Hil. THE NUCLEAR QUADRUPOLE INTERACTION 
Each of the CEF levels may produce a magnetic field and an 
electric field gradient at the nuclear site; this results in hyperfine 
splittings of the electronic levels. A rare earth nucleus thus ex- 
Periences at a certain time a magnetic field and an electric field 
gradient which depends on the electronic state thatis actually popu- 
lated at this time. The situation substantially simplifies at elevated 
temperatures where the spin-lattice relaxation phenomenon produces 
Hapid transition between the different CEF levels. The nucleus 
under these circumstances experiences a magnetic field and an elec- 
tric field gradient which in essence result from averaging these 
fields over all electronic states weighted according to the population 
numbers. This averaging process, which essentially constitutes a 
time averaging process, holds only if the significant electron relax- 
ation times are short compared to all other relevant times such as 
the nuclear lifetimes andthe nuclear precession times, a situation 
prevailing at temperatures above a few degrees Kelvin. In particu- 
lar, the magnetic hyperfine interaction cancels in the absence of an 
external magnetic field and all one is left with is the quadrupole 
hyperfine interaction (23),(24). An example of this situation is 


illustrated in Fig. lb for an assembly of nuclei. The quadrupole 


(23) Rei. Gohen, U.Hauservsand RR. SS eMessbaver coc. 
Mussbauer Conf. 2nd, (John Wiley and Sons, N.Y., 1962) 
li 2 


(24) R.L. Mussbauer, Rev. Mod. Phys. 36, 362 (1964) 





are 


interaction is strongly temperature dependent since the overall CEF 


splitting within the lowest electronic term is only of the order ofa 


few hundred cm > 


| 


The electric field gradient which interacts with the nuclear 


quadrupole moment of a rare earth nucleus bound in an ionic crystal 


has four significant sources: 


1) 


PA) 


3) 


One contribution is the direct field gradient produced at the 
nuclear site by all of the ions surrounding the host ion which 
contains the nucleus in question. This contribution in the 
case of rare earth ions is usually negligible in comparison 
with the contributions from other sources, particularly at 
low temperatures. 

Another contribution results from the electric field gradient 
produced at the nuclear site by the electrons inthe partially 
filled 4f-shell of the host ion. This field gradient results 
from the interaction of the 4f-electrons with the CEF pro- 
duced by the surrounding ions. This interaction efitectively 
induces electric multipole moments (multipole polarization) 
in the 4f-shell; the quadrupole part of this polarization con- 
tributes to the electric tield pradient experienced by the 
nucleus. 

A distortion is usually also induced by the CEF in the closed 
electron shells, yielding another contribution to the total 
field gradient experienced by the nucleus. This contribution 
is proportional to source (1), with proportionality factor “or 


The absolute value of the proportionality factor is in the 








Sita 
case of the rare earths usually large in comparison with 
unity, thereby leading to such an enhanced field gradient 
(antishielding effect) that it often becomes comparable with 
Pre one resulting from source 2). This 1s the “latvice. 
Sternheimer effect (25)-(28). 

4) Another field gradient contribution due to an induced quadru- 
pole moment in the closed electron shells results from the 
miteraction of the closed electron shells with the electroncuam 
the partially filled 4f-shell. This relatively small contri- 
bution, which is proportional to source 2), with proportionali- 
by factor -Ro> is the "atornic'' Sternheimer effect (5), (28). 

Collecting the different contributions, we obtain for any compo- 


nent oot oP tievelectric tield pradient tensor 


Sey. (Lat) 4 _ 
os ) ay (1-Ro) eq. 
where i and Ro are the lattice and atomic Sternheimer factors, 


Mecpectively, as introduced above. 


(25) Peace Wikwer and G, Burns, Phys. Getters 2yec2o ue 
(26) Brice havynge roc: Phys, Soc, 82, 47 (1 Je2) 
(27) R.M. Sternheimer, Phys. Rev. 132, 1637 (1963) 


(28) Pemerecmamand R.b. Watson, Phyo. Nevesioc = 91 0c 
(1963) 
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methe principle axes system of the electric field gradient tenson 


the nuclear quadrupole interaction Hamiltonian Ae associated with 
the (V)-th CEF level of the ion is given by 
(11) 
my, °° 2 a) Ula! (40 _(Lat) 
H Se 1-~R Vv Vv? + 2_y%) + 
Q- 4121-1) \! Ig Ww ie q, Ji3los ue 
(4) a |Y (Lat) (Lat), 11 (72 2 
-Ro) Gg ar - ¢ VD + (= Many ayy Vie (Ly + 12 


mere 1, Ao 1, are the usual nuclear spin oe and Q is the 
nuclear quadrupole moment. The quantities ay nd and (| qt, vy ly 
determine the direct contributions to the electric field gradient at 
the nuclear site produced by the surrounding ions in the lattice and 
by the 4f-electrons of the host ion, respectively. They are defined 
by 


eqit*) = [d2v(r, 0, Y)/ dx, ox, ] (2 


li ree wee 


where V(r, %, Y) is defined in Eq.(1), and 


(13) 


4f-electron 
Lali = 12a? (Feta) [2% 2s) eo] 


where the wave function lv of the y-th CEF level is of the form 
given by Eq. (9). 
Explicitly we obtain for the lattice contribution from Eqs.(}), 


(2), (12) 


2 (Lat) _ 0 2a aay) eet) Z 
. Za eD ae cee Nyy gee ae (14) 


ow 
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Expressing the contributions from the 4f -electrons within a manifold 


of states of constant J in terms of operator equivalents, we obtain 


from Eq. (13) 
a |y> =-Csilels> cry, <v[3s2=24| (15a) 


Cv ay) dy [YD = 3/2) Gel TD Cr ye VISE LE LD 


(15b) 


where Cao is defined by Eq. (7). 

Usually one observes only an average field gradient from the 4f- 
electrons, which is a field gradient from the individual CEF levels 
weighted according to their Boltzman factors, as discussed above. 

If we consider only those electronic states which belong to the lowest 
manifold spanned by the state vector J, then the average direct 
contribution from the 4f-electrons to the electric field gradient acting 


on the nucleus at temperature T is given by (29) 


25+ 1 
f) 
(4f) = <v |g ie Ny» - exp(-E , /kT) 
Sis Pr a (16) 
Sy exp(-E y,/kT) 
v= 1 


The diagonal component of the averaged total electric field gradient 


tensor is according to Eqs.(10) and (16) given by 


(29) A more general description including effects of higher J 
states is given in the Appendix I 
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Le ; 
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where we have neglected any temperature dependence of the lattice 


contribution ee : 

The total Hamiltonian describing the average quadrupole inter- 
action at temperature T may now according to Eqs.(11) and (17) be 
written as 


(18) 


2. ye oo L 2 . 2 
Gee “le (31° ql") - Giese Wy) T 2 (1 ry 


coo 


far) = —<-~ 
! dizi) 


We shall now apply the preceding formalism to the particular 
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ease of Tm 


of n°? (I = 1/2) is not removed by the Hamiltonian, Eq.(18); the 


The twofold degeneracy of the nuclear ground state 


8.4 keV excited state (I = 3/2), on the other hand, is split by the 
nuclear quadrupole interaction into two states. Their energy secpara- 
tion AE... which follows from the diagonalization of the Hamiltonian 
HT) is given by 


“~ 


P2 (19) 


2 : l | 
CAE) p = (e?Q/2) Kaeo + 5 Bee yy TP 


i 
$0 


This expression may be written in more detail, by using Eqs.(5), (8), 


(14),(15), (17): 





See 
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aC 
(OE), =(e?Q/2){ <I a| J? Gs Se rp pe (1-7) 
E 
2 rae SO 
4C ° 
a = f <J|| | i CaO, ae ? or ‘ae Gay | 
1; 


(20) 
where os: = eee and Gs ae oS Der are thermal averages defined 
as those given by Eq.(16), while the parameter cers is defined ay 


Ge ne ee ee 


a a (21) 


> 4f 
It is just this splitting CAE Dp that is measured as a separation of 
gamma lines in recoilless resonance absorption experiments. 

Several additional hyperfine interaction mechanisms which 
Memerabute to the net nuclear quadrupole couplinyp of a rare Garth ion 
have been neglected in our calculations. These additional contri- 
butions arise in second-order perturbation theory with the principal 
effects coming from the magnetic hyperfine interaction itself (30) 
(the so-called pseudo-quadrupole coupling) and from the quadrupole 
interaction with states of higher J.) admixed into the ground state 
multiplet by the CEF. We have made calculations of these contri- 
butions for the compounds covered in this paper and they amount to 
less than 1% of the total quadrupole interaction energy. 

In order to compare experimental results with theory within the 


framework of the CEF model it is convenient to replace in the theo- 


(30) hee eethothwmroc, Pnys, oec. B70), 119 (1957) 





ae 
retical expression for the quadrupole splitting CAE) all quantities 
involving electronic radial integrals (the theoretical determinations 
of which is presently somewhat uncertain) as well as the nuclear 
quadrupole moment by experimentally observable parameters. For 


this purpose we introduce the dimensionless parameters 


p, = eQ <r? CS fla I> /C} (22a) 


P5 Qti-y,)/ GE De (22b) 


Expressed in terms of these parameters the quadrupole splitting 
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in im reduces to 


aa 0 Ce 7 0 aan 
(AED > = | cs ee ee Pa 
(23) 
PENG 2VC. GC oo + dl ean *\2 
3 eee ee an Le 


: : OTe 1S 2 4 2 : . 
The temperature averages < sJU cs ay and oe J* ee within 
the framework of the CEF model depend only on the experimentally 


observable CEF parameters See 





SAE 
IV EXPERIMENTAL TECHNIQUE 

The nuclear quadrupole interaction was measured by using the 
technique of the recoilless nuclear resonance absorption of gamma 
radiation (2). The partial decay scheme of Er Ae is shown in Fig. 2. 
Measurements of the gamma resonance absorption were performed 
as a function of the relative velocity between sources and absorbers. 
The measurements involved sources of erbium trifluoride (Erk) and 


erbium oxide (Er O,) and absorbers of thulium ethyl sulfate 


Zz 
(Tm (C,H,SO,), 9H,O , abbreviated to TmES)and thulium oxide 


(Tm ,O,). 
Anhydrous Erk, provides an excellent source for experiments 
Memazine the 8.4 keV line of Tae The crystal structurecor cme 


heavy rare-earth trifluorides has been investipated by Zalkin and 
Templeton (31). At temperatures below about 900 - 1000° C the 


stable phase is orthorhombic, space group o 


-Pnma, having four 
Memwia units per unit cell. The rare earth ions are crystallogra-— 
phically equivalent, having the point symmetry m. Thus, although 
Meewerectric field gradient tensor ne is not axially symmetric... 
erbium (or thulium) nuclei experience the same oe Therefore, 
the quadrupole splitting of the recoilless absorption line given by 
Eq.(23), may be expected to pass through zero or at least through a 
minimum ata specific temperature (550° K in this case). The advan- 
tages of a single-line source are thereby obtained. The line width 


obtained this way with sources of Erk, is less than with sources of 


(31) A. Zalliam and D-H Penipleton, J. Ame? °Chem, Soc. (>, 
2453 (1953) 
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Er,O, where the presence of non-equivalent erbium sites complicates 
the situation (23),(24). At the same time, with reasonable precau- 


mons, the ErfF, can be maintained at the critical temperature for 


3 
periods of several weeks without decomposition or reaction. This 
chemical stability does not exist with most other erbium salts in 
which the erbium ions are also crystallographically equivalent (e.g. , 
the sulfate, nitrate, chloride). 

Anhydrous Erk, was prepared from erbium metal or erbium 
oxide by a ''wet"’ process. The metal or oxide was first dissolved in 
a small quantity of nitric or hydrochloric acid in polyethylene centri- 
fuge tube. A few ml of aqueous hydrofluoric acid were then added 


and the mixture heated at approximately 100° C in a water bath for 30 


minutes. The somewhat gelatinous Erk, 


precipitate was then centri- 
fuged down, the excess solution decanted off, the precipitate washed 
with distilled water, centrifuged three to five times and dried in air 
at roughly 100°C. Air drying yields a hydrated Erk, of unknown 
composition. To remove the water of hydration, the dry contents of 
the centrifuge tube bottom were transferred to a small tantalum boat 
and annealed in an evacuated fused quartz tube. Experience showed 
that the hydrated Erk, could be converted directly into a mixture ot 
the several forms of oxyfluoride (32) if the annealing temperature 
was raised too rapidly. The procedure finally adopted was to hold 


-5 
the hydrate at room temperature at about 10 © torr for at least 12 


hours in order to pump off most of the water. The temperature was 


(32) W. H. Zachariasen, Acta Cryst. 4, 231 (1951) 
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then raised slowly (in 6 hours) to 150°C, thus removing virtually all 
of the water. Finally, the temperature was raised to 850° C in an- 
other 6 hours and then reduced back to room temperature within 2 
hours. This procedure yielded consistently good clean x-ray powder 
patterns of the orthorhombic phase without a trace of the hexagonal 
phase appearing (31). Erf’, prepared in this manner appears to re- 
main stable at room temperature over an indefinite period of time. 
mueelevated temperatures care must be exgercised to avoid reaction 
with oxygen or water vapor. Sources of Erk, Were prepared iene 
above manner from Er,0, (usually enriched in ae or from erbium 
metal after irradiation in the Materials Testing Reactor, Arco, Idaho. 
Alternatively, the Erk, was prepared first and then irradiated. Iden- 
tical spectra were obtained by the two methods. 

Absorbers of TmF, were used in order to experimentally determine 
the critical temperature at which the narrowest possible emission 
line is obtained with sources of Erk. Figure 3 shows the tempera- 
ture dependence of the quadrupole splitting in TmF,. The source was 
mounted in a small evacuated oven shown in Fig. 4. The absorber 
was maintained ina helium atmosphere within an oven equipped with 
beryllium windows. 

It is interesting to note that the same minimum line width 
(1.8 cm/sec) was obtained in both the trifluoride-trifluoride and tri- 
fluoride-ethylsulfate experiments. This strongly suggests that the 
quadrupole splitting of the trifluoride line does indeed pass very near 


to zero at 550° K (24). This minimum observed line width of 


1.8 cm/sec may be compared with the theoretically predicted line 





Temperature dependence of the quadrupole splitt: 


s 


of the 8.4 keV gamma line of Tn et using an ErF, 


source anda TmF, absorber. Source ana absorber 


were maintained at the same temperature. 
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ae, 4: Details of moveable source oven. The entire oven 
(weight O. 2kg) was moved relative to the absorber 
by the cam drive. The main body of the oven was 
made of stainless steel. A similar oven made ot 
aluminum was used with the transducer drive. The 
heating element was fabricated from nichrome strips, 
1/16 x 0.005 in. For a source temperature of 550° K 


the power input was 30 W. 
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oo 
width. In a transmission-vs-velocity measurement since an emission 
line of width IT is moved over an absorption line of width I, one ex- 


pects a minimum line width of 2f . In the case at hand this is (in 


velocity units) 


et cE = 7 = 0.74 cm/sec 


=o 


based on the lifetime (33),@ = 6.28: 10 Sec, didvenera. 
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ey = 8.42 keV, of the first excited state of Tm However for an 
absorber of finite thickness this width increases by a correction 
factor which is 1.47 in our case for an absorber of 5 mg/cm*% of 
thulium and a total conversion coefficient of 325 from Kankeleit et al. 
(34). Thus the predicted line width is 1.09 cm/sec which must be 
compared with the minimum observed line width of 1.8 cm/sec. The 
observed line width is 1.6 times broader than expected. The origin 
of this line broadening is uncertain. 

Absorbers of TmES were prepared by crushing single crystals. 
Absorbers of Tm,0, and sources of (enriched) Er ,O, were prepared 
from commercially available material. Absorbers of all materials 
to be used below room temperature were prepared by mixing the 


powdered samples with a soft wax and pressing the mixture into thin 


disks between mylar films. Absorbers and sources of all materials 


(33) R.E. McAdams, G.W. Eakins, E.N. Hatch, Phys. Letters 
oy 219 (1963) 


(34) mewisankeleit, fF. Boehm, R. Hager, Phys. Kev. iprescs 








ae 
to be used above room temperature were prepared by settling the 
powdered samples from a slurry of dry acetone onto 1/2 mm thick 
beryllium windows. 

The relative velocities required for Doppler-shifting the gamma 
lines were produced by using both cam drives (35) (providing constant 
velocities) and transducer drives (36) (providing constant acceleration). 
A block diagram showing the experimental apparatus for use with the 
cam drive is shown in Fig. 5. The experimental arrangement used 
with the transducer drive is given in Ref.(36). Proportiona] counters 
filled with one atmosphere of a mixture of 90% argon and 10% methane 
(by volume) and equipped with 1/2 mm thick beryllium windows were 
used as detectors, see Fig. 7. 

Bmeryostat specifically designed for recoilless resonance absorp: 
tion experiments with low energy gamma radiation was used for the 
measurements (37). The sample temperatures in the range from 1O°K 
up to 300° K were attained by either controlled heating of the cooled 
sample holder, by pumping on liquified gases, or by using exchange 
gas cooling. The sample disks were clamped between thin beryllium 
disks soldered to the cryostat sample holder in order to ensure good 
temperature uniformity and stability. Temperature measurements 
were made using carbon resistors and thermocouples. The oven used 
for heating sources to 2000° K is shown in Pieeror 


16 3 
Some typical MUssbauer spectra for Tm ? are shown in Fig, Ge 


5) R.L. Mutssbauer, Proc. MbBssbauer Conf, 2nd, (John Wiley 
and Sons, New York, 1962), p. 38 


(36) E. Kankeleit, Rev. Sci. Instr. 35, 194 (1964) 


(37) F. T. Snively, to be published 








mae. 5: 


Bales 

Block diagram of experimental apparatus for use with 
the cam drive. A detailed schematic of the pro- 
grammer is presented in Fig. 6. All other electronic 
equipment shown is commercially available. The cam 
drive has been described elsewhere (35). The recoil- 
less resonance absorption measurements were per- 
formed by first moving the source at a normalizing 
speed (12 cm/sec), then at a measuring speed (v), and 
finally at the normalizing speed. The three runs were 
of about 5 min each. The counting rates, C, during 
the two normalizing runs were averaged and combined 
with the results of the measuring speed to yield the 
amount of absorption, A(v). 


Ply = | ey sy = yee ae) 


Vs AVE 


This sequence of events was repeated until the counting 
statistics were satisfactory. The information at the end 
of each run was printed out on a typewriter for moni- 
toring purposes and punched out on paper tape tor 
processing by an electronic coniputer as | hesemtiT. 
process was automatic except for changing from one 


measuring speed to another. 
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Schematic of programmer for use with cam drive. 


This programmer controlled the changes from the 
normalizing speed to the measuring speeds and 
initiated the counting and print-out cycles of the 


scaler, “see Tre. 5. 
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Proportional counter. The anode of the counter is 


a 3 mil stainless steel wire kept under light tension 
by a spring in the end fitting. The beryllium window 
is 0.5 mm thick. This counter was operated at about 
2500 V with 1 atm of 90% argon and 10% methane, by 
volume. Under these conditions the resolution at 


8.42 keV was about 17%. 
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-40- 
High temperature, resistance heating, source oven. 
The heating coil was fabricated from two pieces of 
56 mil tungsten wire 15 in long. These two pieces 
were connected in parallel and the source container 
was suspended from the midpoint of each piece. For 
O 
a source temperature of 1970 K the power input was 


4.8 kW. A vacuum of tors torr was maintained in 


the oven. The thermocouple made of tungsten vs. 


tungsten - 26% rhenium is reliable to 3100° K., 
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Quadrupole splitting of the 8.4 keV level of Tm 


in an absorber of thulium ethyl sulfate (5 mg/cm? 


69 in ere 


of thulium). A “single line” source of mia 3 


was used at the critical temperature T = 550° K 
throughout curves a-d. The spectra a, b and c, d 
were obtained by using a cam drive and a transducer 


drive, respectively. 
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The nuclear quadrupole interactions measured as a function of 
temperature in the compounds TmE5S and Tm,03 are given in Figs. 
10 and 11, respectively. Details of the figures are explained below. 
The reduction of our experimental results is carried out in two sub- 
stantially different ways: 

Method 1: We combine our nuclear quadrupole splittings obtained 
from gamma resonance absorption measurements with optically de- 
termined CEF levels and obtain two quantities 

ae (9) t = (c2)"]? and p, (co), + Le 12 which ee 
pend directly onthe electronic shielding factors, compare Eqs.(5), 
(21), (22). This method emphasizes the low temperature data, which 
have the smallest relative errors. 

Method 2: The same two quantities may be obtained without the 
necessity of referring to any optical determination of CEF levels, 
merely by using gamma resonance measurements obtained at elevated 
temperatures. This method is useful in those cases were measure- 
ments can be performed at temperatures which are large compared 
to the overall CEF level splitting, but small compared to the spin- 


orbit splitting. This is the case in both TmES and Tm,0,. 


lL. Thulium Ethyl Sulfate (TmES) 
+ 
All rare earth lattice sites in TmES are occupied by Toe ions 


with point group symmetry C By choosing the proper coordinate 


3h" 
system (22) the relevant CEF parameters as defined by Eq. (8) are 


0 0 


limited to oe Cy: Ce and e. for this symmetry. This leads to an 
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~-45- 
Eemuperdture dependence of the miclear quadrupole 
interaction of ie 2” in absorbers of thulium ethyl 
sulfate (TmES). Sources of ErF, (l= 550° K; single 
line) were used. Curve Ais the best two-parameter - 


fit (parameters Py and p to the experimental data. 


>) 
Curve B is the best one-parameter-fit (parameter P) 
to the experimental data, thus disregarding the lattice 
contribution to the electronic shielding (i.e. o> = ane 0). 
The CEF parameters in set 3 of Table II were used in 
both curves A and B. Observe that (AE) eee) | 
T— oo 

for curve B. 

The difference between curves A and B shows the 
large contribution of the lattice part eo one to 
the electric field cradvemteat the nucleus. Curve 
illustrates in particular, that it is not possible to ob- 


tain a good fit to the experimental data by merely 


adjusting the theoretical value of either Q or Ro: 
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A: fully including shielding 
a a P, = (9.22 + 0.05) X 10°” 
p2* (1.972 0.05) xX 107° 


B: excluding lattice contribution 
to shielding 


P +=6.68 x107® 
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Temperature dependence of the nuclear quadrupole 


169 in Tm,0O,. Sources of Erk, 


(1 = 550° K; single line) and absorbers of Tm,O 
6 2s 


interaction of Tm 


(5 mg/cm?* of thulium) were used for temperatures 

of the absorber in the range ll’ K< T < 700°K. 
Absorbers of TmES (T = 300°K; single line; 5 mg/cm? 
of thulium) and sources of Er,0, were used for 
temperatures of the source inthe range T 7? 700° K. 
Curve A is the best two-parameter-fit (parameters 

Py and p,) to 1thevexperimental data, using the Chr 


parameters of Table VI. The insert shows a typical 


spectrum. The importance of the lattice contribution 


aye) ee to the total field gradient ae ghee 
(1-Ro) a at the nuclear sites 1s strikingly demon- 


strated by the fact, that the quadrupole splitting (AE) 
does not approach zero in the high temperature limit, 
but rather goes through a minimum and then increases 


Ueat) 
again, with ce Ree (I-y) q.; 
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fo 
axially symmetric electric field gradient at the nucleus, i.e. 


cme Boer = 0. In this case the quadrupole splitting <AED,, of 


the gamma lines, Eq.(23), reduces to 


CAE = 20, (0, GL -1D_ + 40,1 (24) 


Method |: In order to obtain the quantities oa and C352 en 
entering in Eq. (24) we use different sets of optically determined 
CEF parameters given in Table Il. Set 1 was obtained for Ta in 
LaES by Wong and Richman (14), who employed observed optical 
levels from a series of different optical multiplets. Set 2 was ob- 
tained for porte in TmES by Gruber and Krupke (15), who again 
used observad optical levels from a series of different optical multi- 
plets. In contrast, set 3 was obtained by a least-squares method 
using only levels observed by Gruber and Krupke (15) in the "Hy 
term of hae in TmES. The evaluation of the c™ given in set 3 
thus does not employ optical terms other than “Hy and therefore 
should be the set most appropriate for our reduction of the nuclear 
quadrupole measurements. To permit a check on the intrinsic con- 
sistency obtainable by using one set of CEF parameters for the whole 
series of optical levels we confront in Table II] observed and cal- 
culated CEF levels. The overall agreement is rather encouraging, 
the average deviations between calculated and observed values being 
only of the order of experimental uncertainties. Table IV gives for 
set 3 the wave functions and field gradients for the CEF levels which 


are necessary for the evaluation of the temperature average 
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CEF parameters ce for thulium ethyl sulfate (units oa 


—— 


Set number References 


Wongp and 
Richman (14) 


Gruber and 
Krupke (3) 


See text 











See 
TABLE HI 


Observed and calculated CEF levels for thulium ethyl sulfate in 
the "H, term of the ground multiplet (units cae The calculated 
levels of set | were taken from Wong and Richman (14). For sets 2 


and 3 the following reduced matrix elements were used (15): 


(sllafs) =1.0197 - 107°, (a |pis) =1.5938- 1077 and 


i] 


(3 ya) = -5.5318~ 107°. 
Observed levels © Calculated levels B 
Seu Set 2 Set 3 
20235 B0G.3 304.7 300.8 
274.0 Zale Lu 2145 
Cees Zoe Die et, 
20229 Zito Zl Se 
r93.9 204. 3 204.0 198.8 
73 eae 161.4 Poe 
O59 oes bile Prose? 
bc) em Zoo 20 
0 -0.5 -4.4 ORF, 
a Optically determined levels of Gruber and Krupe (15) 
b Calculated levels using the CEF parameters given in Table UL. 


The center of gravity of the calculated levels is adjusted to 
give the best fit. 
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TABLE IV 


Energies, wave functions and electric field gradients of the CEF 
levels of the oH, term of the ground multiplet of thulium ethyl] sulfate 
(C,, symmetry), using the CEF parameters of set 3 and the reduced 


matrix elements given inthe caption of Table II. 


Energy Degeneracy Wave Function * 35% al 
7. 1 -0.707|-3) + 0.707 |+ 3) -15, 
110.6 2 -0.446|-2) + 0.895 |+ 4) a 


0.895|-4) - 0.446 |+ 2) 


58.0 1 0.697 |-6) -0.168|0)+ 0.697 [+ 6) 63. 
51.8 ! -0.707|-6) + 0.707 |+ 6) 66. 
35.1 2 -0.305|-1) + 0.953 |+ 5) 26. 


-0.953 |-5) + 0.305 |+ 1) 

5.9 ! Ores) ae GO. 7Orr ts) -15, 

73: 0 2 0.895 |-2) + 0.446 |+ 4) =a 
0.446 |-4) + 0.895 |+ 2) 

131.7 2 0.305 |-5) + 0.953 |+ 1) yD, 
0.953 |-1) + 0.305 |+ 5) 


-163.0 ! 0.119 |-6) + 0.986 jo) +0.119{+ 6) -39. 


a The general form of the wavefunction is given in Eq. (9) 


ae 





= 
C352 - J*>,., Eqs. (16) and (24). 

A summary of the reduced data obtained by combining our 
measurements of the temperature dependence of the nuclear quadru- 
pole interaction for ere in TmES with the results of optical 
measurements performed on the same compound is presented in 
Table V. The dimensiongless parameters Py and P, presented in 
Table V are experimentally obtained quantities (compare Fig. 10 
curve A) which hold within the framework of the CEF model. The ad- 
vantage of introducing these parameters is that their deduction does 
not depend on a knowledge of the radial distribution of the 4f-elec- 
trons or the value of the nuclear quadrupole moment. Such a know- 
ledge, however, enters into the evaluation of the shielding factors, 
fas. (22), (21) and (5). 

Method 2: At elevated temperatures the temperature average 

Gi? = Je entering in Eq.(24) may be approximated by Eq. (I-17) of 
Appendix I, which yields (352-5, =-14.1 CO/KT. Expansion 
(I-17), which applies to the case of an axially symmetric field gra- 
dient, was first given by Elliott (39). Details are given in Appendix I, 
which also includes an extension to the case of non-axially symmetric 
field gradients. . From a plot as a function of 1/T of our measure- 


ments obtained for TmES at temperatures T > 200° K, Big. [2 we 


é ~ 
obtain from Eqs.(24) and (I-17) the values p (Co) =(0, 182 0/02)cu ‘ 
and Poe = (2.8+1.1) None ee These values may be compared 


(38) M.C. Oleson and B. Elbek, Nuclear Phys. 15, 134 (1960) 


(39) R.J. Elliott, Rev. Mod. Phys. 36, 385 (1964) and private 
communication 





(g¢) Yeqiq pue uoss[O Wor Eg[-Wy jo 
22eYS AFH FS FY JO JUoWIOW ajodnaipenb ieaponu 9yj 103 UTeG G*T = © 9NTeA [eoITJaIOVYY 94 SsuIsy eke 


nnn 


061 on -.01 + (20 ¥ 01) mer 6 0F 2i)- 2 *oruy 
L8} eo! -(800 OF 110°1) ¢- 01 -(€00 ‘OF 191) — : I Poe Ul, 
06€ el ¢-0! woes urs 1° Z) 9-01 -(¢€# Ol) € 2 SAUL 
89€ (oon gO - (S0°O # 26°T) 9-0! . (S0°O = 272 °6) € I SqGwyL 
SOF ZO os (SO°O = LIZ) g-Ol - (SOO ¥ L0°6) 2 l SAUL 
88¢ 1°01 mee (SOO * 80°2) g-01 . (S00 F PE 6) I 1 SquL 

I 

an ae). / Sake 1) 9-7) eg ly aS pouleyFw OC punodwi04y 


I 


‘sqUaUaInseau [E91IdO |Y3 UT SATJUTEJIVOUN apNy{ouT you Op pue UOTE S{UY 

ajodnipenb i1eafonu ay} jo syuauaInseau Ino Wodrj Zutstie sioOiia ATuO aie ¢ pue p SUUNTOD UT pees 
SIOII9 sUy *aTqey Fu} NOYSsNOoAY. pesn aie sjtun 9IWOWW “el ‘81.q areduiod ‘AzIptTea Atey, Jo esuel ayy 
ul ] xtpueddy jo suoteurtxoidde ayy sutsn ‘squiod ejep [equawtiadxa ayy 0} (¢7) OF Buiy3yty Aq pauteiqo 
azam °d pue Ig saajaureaed ayy 7 poyjyaw jo ased ayy ut ‘ATieTMWIg “[I pue Ol ‘s81.q ul uaats squtod 
eyep jequournazadxa ay} 03 ‘(¢2) bg “Way) JO 41j sarenbs-yseat e Aq | poyjyeu jo aseo ay} Ut paute qo 


aizam ¢d pue Tg szojauresed SsaTUOTSUaUIIp 94 J, tofu uI pue Sqwy ut Ulf, IOJ eyep pa onpsay 


tae 


A ATaEVL 





(fe F 
iy by 
‘em -_ 


169 


Nuclear quadrupole interaction of Tm 


| in absorbers 
of thulium ethyl sulfate and thulium oxide plotted as a 
function of 1/T in the high temperature ranges where 
method 2 is applicable (see text). The straight lines 


are the best fit to the experimental data points. 
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Ge 


with the values obtained by method I: Pre: = On Moe O Only) eee 


: =2r 1 = 0. 0) 107? aoe The agreement between these 


and p40 
two values obtained by two different methods gives confidence in the 
consistency of our analysis. In particular we conclude on this basis, 
that our results are not seriously influenced by any temperature de- 
pendence of the CEF parameters eae within the temperature range 


studied (9. 6° K - 300° K). The agreement obtained for the results of 


methods 1 and 2 indicates that the CEF parameter C) is reasonably 


0 
6 


and eo should be even less dependent on temperature, because of the 


independent of temperature. The higher order parameters oa C 


faster convergence of the associated lattice sums. The justification 
of the neglect of any temperature dependence of the CEF parameters 
c™ in our analysis is supported by measurements of Gruber and 
Conway (40), who determined by optical methods the energies of CEF 
or) _ O Ca- O 
levels of Tm tons in PimEo at T - 7% K, 194 K and 273K.” the 
changes with temperature in the position of the levels typically are 
less than 10 em. We therefore feel justified in using in our analysis 
one set of CEF parameters, Cc”, determined optically at a single 
temperature. 
2. yehnuwlium Oxide (Tm ,0,) 
Sr 3 . 7 
see Im ions in the Tm ,O0, (space group T,) occupy two non- 


equivalent lattice sites; sites with symmetry C, and C3. occur in 


(40) J.B. Gruber and J.G. Conway, J. Chem. Phys. 32, 1178 
(1960) 
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the ratio 3:1. Experimentally, only the higher populated ionic sites 
associated with point group symmetry C, are observed. By choosing 


a proper coordinate system the relevant crystal field parameters for 


are 0 2 0 2 -2 4 -4 0 
C, symmetry are limited to Cs C,. Cy Cy Cy : Cy & , Ce: 
er, me, ey ; oe on Gat The parameter a has been elimi- 


nated by the proper choice of the coordinate system. One may always 
a -m . m 
eliminate one Ch by such a choice (22) when C. also exists. 
. : m -m | 
Assume that in a given coordinate system both Cc. and Ch exXUS6, 
Then by rotating the coordinate systems about the Z-axis by an angle 


6 defined by 


one may reduce the problem to one parameter in the new coordinate 


’ 
system, oe which is related to the old parameters by 


No eae: 


n 

The quadrupole splitting of the gamma lines produced by the non- 
axially symmetric field gradient is given by Eq. (23). | 

Method 1: Gruber et al. (16) have studied the optical absorption and 
emission spectra of ea ions in Y,0, at the C, Syinmietry sites, 
the same CEF levels were obtained in preliminary studies of Tae 
in Tm 03, within the limits of the experimental accuracy. Using 
the energy levels obtained for diluted eae by Gruber etal. (16) we 


have calculated the crystal field parameters c™ which are included 


in Table VI. This calculation was performed by minimizing Xs 
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TABLE VI 


+ 
Observed and calculated crystal field levels for Teena in thulium 
oxide at sites with C, symmetry, inthe “Hy term of the ground 


multiplet (units enh The following set of CEF parameters was 


0 | ae 0” ae a 
Ae ware 0 | oe ee 4 
a ae 6 | -6 


Calculated Observed Calculated Field Gradients 


A preliminary set from Gruber et al. (16) 


From Gruber et al. (16) 


eves cei: eo ey (3? — J?) Cisne) 
770 796.9 ar -74, 
768 7eaee - 4, mar 
680 ee Aan 
674 iP ae 
497 494.0 ! eae 7p 
429 Aa5 7 ! oe 17. 
344 ae ie. 
286 340. 0 l a 25, 
258 230.3 ! =e) fie 
200 219.0 ] ce oe 
95 89.3 l Lao) ae 
44 30.7 ! 40. 28. 
at 0 ee 37, 








FEO 

using the technique described by Davidon (41). A summary of the re- 
duced data obtained by combining our gamma resonance measurements 
for tha in Tm,9, (compare Fig. 11) with the results of the optical 
measurements is included in Table V. 

Method 2: At elevated temperatures the temperature averages 
<1--J*>., and <ag + Le entering in Eq.(23) may in first order 
be approximated by the expressions given in Eq. (1-17) and (I-18) of 
Appendix I, which yield Gia i =-14. 1] Co/kT and Gir i py 
=o. 4 CS/KT. From a plot as a function of |/T of our measurements 
obtained from Tm,90, at temperatures T = 127.0 Ee Fig. 12, we ob= 
tain by using Eqs.(23), (I-17), (1-18) the following values: 

ac. [ic%) ' ey }2 = (0.5£0.3)cm~2 and p, ((C%) + Licey = 
(4.0+0. 9) ie aes cm7} These two values again may be compared with 
the corresponding values obtained by method 1, namely (0.50840. 001) 
ae ema ( 3, 80 + 0.03) Nore cm” !, respectively. The rather good 
agreement between the values obtained by methods | and 2 suggests, 
as inthe case of TmES, that the neglect in our analysis of any 
temperature dependence of the CEF parameters co is a justifiable 
approximation. Furthermore, the agreement suggests the absence of 
crystallographic phase transitions in the whole temperature range 
studied. X-ray diffraction studies of Stecura and Campbell (42) do 


not reveal any phase transitions within the temperature range 300° K 


<T< 1568°K. 


(41) W.C. Davidon, ANL-Report 5990 rev., (Nov. 1959) 


(42) S, Stecura and W.J. Campbell, U.S. Bur. Mines, Rept. on 
Investigation No. 5847 (1961) 
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Vir BPEEGPRONIGC SHIELDING PACTORS 
Our experiments reveal the presence of strong charge polari- 
zations of closed electron shells. The shielding (or antishielding) 


factors R and o>, which were introduced in Sections II and III, 


Q' Yo 
are a measure of these charge polarizations. 

The antishielding factor a ("lattice'’ Sternheimer factor) may be 
calculated by several techniques when the free ion wavefunctions are 
known. Wikner and Burns (25), Ray (26), Sternheimer (27), and 
Freeman and Watson (28) have made calculations of this quantity for 
certain rare earth ions, and their results are summarized in Table 
VII. Wikner and Burns used the (restricted) Hartree-Fock wave- 


- 
functions calculated by Ridley (43) for oye! and ea and calculated 


y.. by means of a perturbation-variation method. Sternheimer used 


QD 


the same wavefunctions, but calculated i by direct solution of the 
inhomogeneous SchrUdinger equation for the perturbed wavefunctions. 
Freeman and Watson used the unrestricted Hartree-Fock formalism 
to calculate es for Cote. The value of Freeman and Watson for 
ee is not very different from that obtained by Sternheimer for the 
neighboring ion ee , but differs appreciably from the value which 
Wikner and Burns obtained for ar 

Theoretical evaluations of the shielding factor Rg are Imnore 
involved. This results because of the proximity of the closed elec- 
tron shells to the distorting source, the 4f-electrons. For this 


reason, one may even expect that the distorted shells may produce 


(43) E.C., Ridley, Proc. Cambridge Phil. Soc. 56, 41 (1960) 





269. 

repercussions upon the 4f-electron shell, as was pointed out by 
Freeman and Watson (12), (44). Table VII includes the few available 
theoretical values of Ro = ead + Dane for rare earth ions. 


Theoretical evaluations of the shielding factor o0, are physically 


L 


similar to those for Ve The additional complication arises from 
the fact, that ce is a measure of the closed shell distortions ex- 


perienced at the origin, while o, is a measure of the closed shell 


2 


distortions experienced at the position of the 4f-electrons, thus re- 
quiring in addition a rather precise knowledge of the 4f-electron 


density. Theoretical predictions for 0, are still rather qualitative. 


ie 


Lenander and Wong (10) came to the conclusion that the shielding 


factor o 5 was of the order of 0.5 to 0.75 in the case of PrCl, while 


Watson and Freeman (12) in the case of cerium ions likewise con- 


cluded that shielding via the o, factor is large. Ray (26) arrives at 


Z 
the theoretical value of c, =O. Demorutme case Of PrCl,. Burns (9), 
on the other hand, using analytic perturbation calculations, con- 
cludes that the shielding factor o5 should be at most of the order of 
Oeetor rare earth ions. 
The analysis of our experimental results yieldsthe parameters 

Py and P, given in Table V. Using the value of the nuclear quadru- 
pole moment Q we can evaluate the parameters — and (lay eae 
>... Values of these parameters are included in Table V. It 


appeared reasonable to use a theoretical value for ve to obtain the 


radial integral Cr: since theoretical evaluations of this quantity 


(44) A.J. Freeman and R.E. Watson, Phys. Rev. 131, 2566 (1963) 
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TABLE VII 


— 


Theoretical values of Sternheimer shielding factors for rare 


earth ions 


Ion Y R 


oO rad 
Ge = 735 Owe 
Ee - 16.4 
ee” -105 
e - 78.5 
eA - 61.5 
7a -~ 74 
Bae 


a Using the value (28) (r-*) 


4f 


ang 


=.4,7l a.u. 


Reference 


(28) 
(26) 
(25) 
(27) 
(25) 
(27) 


( 5) 
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mppear to be relatively reliable. 

Besides the values of the "electric'' radial integrals < © and 
ee which enter in the analysis of our measurements of the 
nuclear quadrupole interaction, there exists a "magnetic" radial 
integral <r, which enters inthe analysis of nuclear magnetic 
interactions. The effective integral (ry, likewise may be asso- 
ciated with a shielding factor (5), (45), which in analogy with the 
electric case is defined through the relation (ry, = Gare. 
(compare Eq.(21)). The difference between the values of a. and 
©), arises because the contributions from the closed shells differ 
for the quadrupole and the magnetic interactions. This difference is 
due to the different forms of the interaction operators for the nuclear 
quadrupole, magnetic orbital and magnetic spin interactions, as was 
emphasized by Sternheimer (5), (45) and Freeman and Watson (44). 

The radial integrals Ge Cae and cD oe which enter the 
nuclear quadrupole, nuclear magnetic and CEF interactions, incor- 
porate the contributions to these interactions from both the partially 
filled (4f) and the closed electron shells. These radial integrals in 
principle may be taken from experimental observations, a procedure 


adopted inthis paper. Table VIII includes a compilation of relevant 


(45) Re ec men nem nye. (eV eOr a5) OM a2) 

(46) im Lindgren, Nuclear Phys. 32, 151 (1962) 

(47) Pee deanna Io ulanderen, Phys. Kev 22, 1802 (1961) 
(48) Ree ceonen, E hyc.wkev. 134., A94 (1964) 


(49) E. Gerdau, W. Krull, L. Mayer, J. Braunsfurth, J. 
Heisenberg, P. Steiner, E. Bodenstedt, Z. Physik 174, 389(1963) 





Shige 

Sites ; . ; 
radial integrals for Tm ions in different chemical surroundings, 
obtained from experimental data by using the theoretical values for 
Q and De given in the caption. The table includes radial integrals 
evaluated from our gamma resonance studies as well as from other 
pertinent experiments. 

The interpretation of the radial integrals interms of electronic 
shielding factors requires a knowledge of the quantities wane and 
(2?) as discussed above. These radial integrals are not acces- 
sible to direct experimental observation and one is forced to use 
theoretical values, the evaluation of which is presently somewhat 
uncertain because of the lack of sufficiently accurate atomic wave 
functions for rare earth ions. Any evaluation of electronic shielding 
factors is therefore limited by the uncertainties in these theoretical 
values. Nevertheless, by using a specific set of theoretical values 
throughout the whole analysis, one still can observe the general trend 
in electronic shielding. 

Freeman and Watson (50) discuss the theoretical situation in the 
evaluation of (x7?) and (r?) for most rare earth ions. These 

4f 4f 
: . -3 

authors, in particular, have shown that the values of Cr ie for 
rare earth ions incorporated ina solid do not differ very much from 
the free ion values (28). Table VII includes a compilation of elec- 

: ae: . 
tronic shielding factors for Tm ions obtained by using the theore- 
tical quantities given in the caption. The uncertainties of the theore- 


tical values of Q, Cai. and (ry, are presumably less than 30%. 


(50) A.J. Freeman and R.E. Watson, Phys. Rev. 127, 2058 (1962) 
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It appears from Table VIII that the shielding of the nuclear qua- 
drupole interaction and the nuclear magnetic interaction, expressed 
through the shielding factors Ro and Ry is always small in the case 
of tm ions. 

As concerns shielding factors other than R, we note again that 
our experiments provide only the ratio Cha yes) Coe = 
(Imy J/{O-«,) Cra ], compare Eqs.(22), (23). It appears from 
column 7 of Table VIII that there is a substantial electronic shielding 
associated with the shielding factor o>, which describes the fact that 
the 4f electrons do not interact with the direct CEF, but with a CEF 
shielded by core electrons (primarily Soa electrons). This ex- 
perimental abservation is in qualitative agreement with conclusions 
we draw for praseodymium salts from NMR measurements on lantha- 
num salts ee ee by Edmonds (6). CEF shielding effects of 
comparable magnitude were also obtained by Blok and Shirley (8), 
in the case of several rare earth ethyl sulfates and rare earth double 
nitrates, using nuclear alignment techniques. We note that our con- 


clusions concerning o, are in agreement with the theoretical 


Z 
estimates of Lenander and Wong (10), Ray (11) and Watson and 


Freeman (12), but are in serious disagreement with theoretical 
conclusions of Burns (9). 


Hels interesting to note the difierence in the ¢_. values presented 


2 
in Table VIII] for TmES and Tm,O,. This seems to indicate that o>, 
depends on the chemical environment, which might result from 


different amounts of overlap of ligand wavefunctions with the central 


rare earthion. This seems to conform with similar conclusions by 





ae 


Hutchings and Ray (21) for the case of PrCl, and PrBr,. It should 


be observed that our conclusions concerning g, are based on the 
plausible assumption that ae is much less dependent on the chemical 
bond fan o>: 

Furthermore, we emphasize that we have neglected the 
non-linear shielding effects (12) in our analysis. Appreciable non- 
linear shielding would invalidate the CEF parameterization scheme. 
However, due to the overall agreement reached in our analysis -- in 
terms of linear shielding -- of the optical data and our quadrupole 
data we conclude that non-linear shielding effects play only a minor 
Pole. 

Similar measurements of the nuclear quadrupole interaction in 
TmES as those reported in this paper were performed by Hufner 
etal. (51). HUfner et al. in the analysis of their data did not take 
into account, that the optically determined CEF parameter C, docs 
not represent an unshielded CEF parameter but rather represents a 
potential at the 4f-electron sites which undergoes shielding by closed 
electron shells of the order of 70%, as shown inthis paper. We 
would like to emphasize, in this context, that the lattice contribution 
to the total electric field gradient at the nuclear sites is most easily 
Soaerved in the measurements performed at higher temperatures. 

The nuclear quadrupole interaction in Tm,O, has been investi- 


Za 


gated previously in a limited temperature range by Cohen et al. (23) 


(51) S. Hufner, M. Kalvius, P. Kienle, W. Wiedemann, 
PeePichienry) 2. ft aysik. 1/5, 416 (1963) 





ae 


and Kalvius et al. (52) using recoilless resonance absorption of 
gamma rays. Inthe analyses of these papers the direct contribution 
meerm the lattice to the electric field gradient, enhanced by electronic 
shielding, was not considered. Preliminary results for Tm ,0, ina 
limited temperature range were reported by Cohen (53). 

Although the importance of shielding effects expressed by the 
factor o, 1s well established, the absolute values of the shielding 


Z 
factor 7, may be in error by up to 30%. On these grounds we do not 
feel that there exists any of the serious discrepancies reported by 
Hufner et al. (51) between the value of the nuclear quadrupole mo- 


ments obtained by gamma resonance measurements and those derived 


from Coulomb excitation measurements. 





(52) Mr Kalvius, W. Wiedemann, R. Koch, FP. Kienle and H. Eicher 
Za Pavysvic 70, 9267401962) 


> 3) R.L. Cohen, Ph.D. thesis, Department of Physics, California 
Institute of Technology, Pasadena, 1962, unpublished 
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VII. SUMMARY 

This work has demonstrated that the technique of gamma reso- 
nance absorption provides a sensitive method for investigating elec- 
tronic shielding by closed electron shells in rare earths, via measure- 
ments of the temperature dependence of the nuclear hyperfine inter - 
actions. It was shown, in particular, that in those cases where 
measurements can be performed at elevated temperatures one can 
obtain information on electronic shielding factors without the necessi- 
ty of relying on CEF parameters determined by other methods, such 
as optical spectroscopy. Our results lead to the conclusion that the 
distortions induced in the closed electron shells by the 4f shell only 
produce a small shielding of the 4f electron contribution to the total 
field gradient at the nuclear site ("atomic" Sternheimer shielding 
factor IRo| < 0.1). Onthe other hand the distortions induced in 
the closed electron shells by the CEF lead to substantial enhance - 
ment of the direct electric field gradient produced by the surrounding 
ions at the nuclear site ("lattice'' Sternheimer antishielding factor 
ies) as well as toa substantial reduction of the CEF as seen by the 
4f-electrons (shielding factor o,). We obtain values for eee 
(l-o,) Otc 590) fOr ee ions in thulium ethyl sulfate and of 128 for 
ail in thulium oxide. The difference in these two values seems 
to demonstrate a dependence on the chemical bond. 

It is interesting to note that the ratio of l-o, for TmES to l-o, 
for Tm,9, agrees approximately with the ratio of the overall] CEF 
splittings in these two compounds. 


It appears that measurements of the nuclear quadrupole inter - 





are presently much better suited to determine electronic 


| ling factors than to determine nuclear quadrupole moments, due 


e existing uncertainties in the different electronic shielding 








See 
APPENDIX G1 
In what follows we derive an approximation for the electric field 
gradient which holds at elevated temperatures. 
The relevant matrix elements entering the expression for the 
electric field gradient introduced in the text, Eq.(15) can be ex- 


pressed in terms of spherical harmonics by 


(3|Jo|| 3) (352-5), = av TS (= Yoder) = (I- 1a) 


(I- 1b) 
(5 lols) (a2 + 52). = NB IB (EL YS (9., 0.)+ ¥5°( 9. €.)]) 
+t -/ T : Canara c Pot si 
where the z extends over all 4f electrons. 
Using the density matrix formalism the thermal average of 


the spherical harmonics in Eq. (I-1) may be written as 


<<) | ee 2, (AM |Y" exp [- B (Hy + V)}|AM) | (ian 
where 
Z = Ss (AM | exp 3 a V)} {aM) (I-3) 
1M 


and B = 1/kT. The Coulomb and spin-orbit interactions are 


represented by Hy and V is the CEF potential defined by 





vas Ga Aeon, (1-4) 


and the electronic wave functions |r v) for pure Russell-Saunders 


coupling are defined by 
(Ho + Vv) [av) = (E, + Ey) jay) (1-5) 


The quantum numbers a@LSJ are represented by A and V isa 
quantum number characterizing the CEF levels. Since the trace of 
a matrix is invariant to the choice of basis functions, we choose 
eigenfunctions of es in Eqs. (1-2) and (I-3) rather than using the 
eigenfunctions in Eq.(I-5) which are mixed in M. The "lattice 
sums'' Qa" introduced in Eq. (1-4) are linear functions of those 


used in the text (compare Eq.(1)). We have for instance 


A> =4Vr TS Ay (I-6a) 


2 -2 Z 
4 -_ 
Bie Oss 4.) 24-715 A, (I-6b) 


and we choose Ao = 0. 
According to Van Hove et al. (54) the exponential factor in 


Eqs.(I-2) and (I-3) may be expanded as follows 


(54) L. Van Hove, N.M. Hugenholtz, L.P. Howland, Quantum 


Theory of Many Particle Systems (W.A. Benjamin, Inc., 
New York 1561) p. 82 
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n= VU 


p B, Boy 
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-exp [ eee e5) tt] V...V exp{ -( Bete ee Vexp aie Hy) 


For a temperature large compared with the CEF interaction energy 
fia. e. pi, < 1) only the first few terms of Bq (1-7) need be considered: 


Hence Eq. (1-2) reduces to 
(1-8) 
—)- eee de CAM | ioe | x'M!) CA'M! ee p, + O(p") 1AM) 
A 
MM! 


where 





iy M O (p* 
24a Jp te, + O18") |AM) 


and 








a7 Ge 


QM! a |AM) = exp (- BE} 6 yy: VEE ; 


- Bexp(-BE ,) (AM' Iv|aAM) for a=at 


»'M! 1M) = 
| exp[B(E .,-E,)}-1 


E.-E 


~exp(-BE (aim! |v |aM) 


a) 
r 


for x + 


Furthermore, if the temperature is also small compared with the 


represeners 


Spin-orbit splitting (i.e. B( EY -E_)>>1, where hy 


1 Q 


Mvemtinst excited term of the ground multiplet and 4A. represents 


0 
the ground term) only the ground term is appreciably populated. 


Thus Eq. (I-8) reduces to 


at = -27) Dy 


De |e agM [YE [ag M) (Ag M'|V[ayM) + 


(i= 


+2 aga [ym |arme) Grae |v age) (ey ] 


' 
M' FXG 
Hoon >O,. Here Z = 2J yt 1 and we have chosen 7 =). Further 
0 
more because of the Wigner-Eckart theorem and the properties of 
the vector coupling coefficients (the notation of Edmonds’ (55) is 


used) we have used the following relations 


(55) A.R. Edmonds, Angular Momentum in Quantum Mechanics 
(Princeton University Press, Princeton, 1957) 
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ze (aM|y"|,M) = (an+ 1)7? (x fy fad> = (-1)" "My ma-M |sanm) 
M M 


) 


(ete See 1)? (IMJ -M|JJ 00) 


z (33 00| JMJ-M) (IMJ-M|JJ nm) =e eee 


and therefore 
z CAM | ¥o"|AM) = 0 for n> 0 (121 


According to Eq. (I-1) we are only interested in the cases of 


even n and im, for which we obtain from Eq. (1-9) 
(9.9 Iq =H (21,t AL™ { Grol, (9) [29) - 


ZAC) (ro 169) [ro B+ 2 ~ (a9 LY 9) [aD 


t 
Xr Ag 


\— | 
paar. )¢ x Jr J.) | Ao) (E,,)7']}- (Qm+i) (i121 


In arriving at Eq. (I-11) we used Eq. (I-4) and the following 


relations: 





Bays 


2a (doM| YE | aM") (atm! | v3 |agM) = 


= [(2n+1) (apt 1)]? (ro t¥, fad a IY, | Ag) 


meee (J) J'nm|J,MJ'-M') (J'M'J.-~M|J'S, pq) E 


= [(2n+ 1) (2p+ 1)] “2 (y, Jy flay Ca IY, apt ye 


6 6 A(JgJ n) 


np m-q 


where A(J,J'n)= 1 if Jo J' and n satisfy the triangular condition 
and A is zero otherwise. 


Following Elliott and Stevens (18) we now make the following 


correspondence between reduced matrix clements 
(Ila ls) =4V4 75 (aLss [> ¥, Y,( 9.) Jalss) Na, (I-12a) 


] 


(la a+ t= -aV*75(aLSS |» Y,(o)feLss +1) INR, (12d) 
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(I-l2c) 
where 
evs erel). (2) 4. lee oat +: 3) (I-13a) 


yg 





ee 


Q =3 J (J+ 1) (st Die see yee et 3) (1-13b) 


a tl 


ze 
i Ci VL) ee en et) 4 to) (I-13c) 


Finally, by combining Eqs. (I-1), (1-6), (I-11) and (I-12) we obtain 


the following expressions: 


(Sfe|s) 333-2 )q 


AY Ge Bag) (I-14) 


I] 


(SJels) Gh +L)» 


AS (r2\_ O(T) (I-15) 


where 


-] 5 fio ffs) |? aupiee 2|(3 fo s*1) [295 54, 


kT ee] 


O(T) = -z (23 +1) 


2 a ffe fa 7) | On am 


' Se (1-16) 


) : 
The factors re] and Qs 52 


from Eqs. (I-13b) and (I-13c) by changing J to J-1 and J-2, 


in Eq. (I-16) are obtained 


respectively. The energies v4 and Ey $2 are those of the 


center of gravity of the terms nearest the ground term which have 


Guantum numbers J+] and J+2. 
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Eipplying these results tothe case of Tm  , the effect of the 


second and third terms of Eq. (I-16) is negligible ( < 1%) at all 


temperatures used in our experiments. Under these circumstances 


we arrive at the following high temperature approximations used in 


i 


method 2 of the text (compare Eqs. (23) and (24)) 


~ = (23+ 1)7' co (alos) ae [xT 
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APPEND it 

Most of the computations involved in the analysis of the experi- 
mental data presented in this paper were carried out with the aid of 
an IBM 7094 computer. This section contains the listings of the 
Fortran IV computer programs that were used. The analysis involved 
four major steps. 
1) OPTIC ... A least-squares fitting of the calculated CEF levels, 
using the method of Section II, to the optically observed levels in terms 
of the CEF parameters. 
2) ICARME ... A calculation of the reduced matrix clements (compare 
Eqs.(4) and (I-12) in the intermediate coupling approximation starting 
with the relevant Slater integrals and the spin-orbit coupling para- 
mreter. 
3) QTAVE ... A least-squares fitting of the calculated nuclear quadru- 
pole splitting, using method | of Section V, to the observed splitting in 
terms of the parameters P) and p>. 
4) HITEMP... A least-squares fitting of the calculated nuclear quadru- 
pole splitting, using method 2 of Section V, to the observed splitting at 
migh temperatures. 

Some of the programs used are available in the SHARE library 


and are not included here. 





ORG 


The required card decks are: 


MIN 
READY 
AIM see SHARE no. 980 
FIRE ZO ANFZO13 
DRESS 
Se 
MATMPY 
PCN 
THEORY 4 
QSQUAR 
ll. LEIGEN 
ie HE RM see SHARE no. 884 PK HMEE 


13. RDM see DEAKE No. 1359 G5 AGC Weg. 


Decks 1-7 were converted from Fortran JI to Fortran IV and 





deck 12 was modified so that it could be used with a Fortran program. 


The listings of decks 8-12 follow. 












_— we ee ee ~~ em eee ee . ~ 


_ ee ee = 


me oN ... 
‘THIS SUBROUTINE I[S THE LINK BETWEEN THEURY4 ANDO QSQUAR TO MIN, 
THE VARIABLE METRIC MINIMAZATION ROUTINE, SHARE NO. 980 
— INPUT DATA | oe _ . . a 
PAR(T)=CEF PARAMETERS IN THE FULLOWING ORDER C20,C22,640,C42;, 
C4-24044,04-4,060 C62 4C6- 24064, C6-44C665C6-6 OR C20,C4U,C6U, 
C66 
— NE=NO. OF ENERGY LEVELS MEASURED 
ITIME= MAX. XEQ TIME ALLOWED 
ALOWER=LOWER LIMIT OF RELATIVE OLFFERENCE USED IN CALCULATING 
GRADIENT 
UPPER=UPPER LIMIT OF RELATIVE DIFFERENCE USED IN CALCULATING 
GRADIENT 
E(I)=MEASURED ENERGY LEVELS DF 8 = 
R(T )})=ERROR [N MEASURED ENERGY LEVELS 
VOUTI=LEVEC-1DENTIFICATION = rr 
OP{IT)=STEP SIZE USED FUR GRADIENT 
| AB=IDENTIFICATION FOR PARAMETERS 970000000 TTT TTT 
SUBROUTINE FCN(NP,G,F,PAR,MFLAG) 
COMMON/COMTHE/NFORM,P(3,40) — CS eee 
COMMON/COMQSQ/NE,E(100),S(120),R(100) 
SecaMMON/COMMIN/H 
' DIMENSTON H(40+40)46(40) ,PAR(40),V(100),0P(40),AB(40) 
Be DIMENSION AH(50.22)48H(50,1) FP (40),1TER(40) pL TERM(40) »DEV(100) 
[F(MFLAG-l)lels2 
“1 READ(5+100)} NE, ITIMEs,ALOWER, UPPER 77 ae 
100 FORMAT(I3 /16/2E12.5) 
Meecret ICLOCK ({iNTIMEey 
[TIME=ITIME+INTIME 
fee SOO)LE(TIGR(T)yV{ by el=leNE)} a 
800 FORMAT (2F20.5,A6) 
_ READ(5,200)(DP(1), IT=1,NP) he a = © an = 
200 FORMAT(F20.5) 
es, 500) AS © = = ae a 
500 FORMAT (1246) 
READ(5,900) 
900 FORMAT(80H HEADING 
_ l ) 
WRITE( 6,700) CO Palle ji, I=1 NP) 
— 100 FORMAT(39HOSTEP SIZE USED IN CALCULATING GRADIENT/(3HO 8€14.5)) ~ 
WRITE(6,600) AB 
600 FORMAT(23HOORDER OF PARAMETERS X/(20A6)) 
NODEG=NE-NP 
~ 2 NFORM=O x | — — ; 
IF(MFLAG.NE.1) NFORM=]1 ; 
LRG) Cs SS) a 
3 P(2,1)=PAR(1) 
— -F20.5#QSQUAR(1,2,2) 
NFORM=} 
DO 10 [=],NP 
me __ITERM(1)=0 
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ITER(I)=0 

Mel 2 P(3,1)=PAR(IT)#(1.0+40P(1)) | _— ea ee — 
FP(T)=O0.5#*QSQUAR (15312) 

Meee ABS({(FPCIL)J-FISTPPAIY+F)) 
[F(RF.GT. ALOWER) GO TO ll | 


mM  4OP(I)= 5.0#DP(T1) a : 5 ee ae : ial 
ITER(IT)=ITER(I) +1 
—_— MeeeereR(1).L1T.93..6G.10 12."  —." © —— 
FP(I)=F 
Meer ieé (6,2001) T,DP OCT) 0 
2001 FORMAT (35HOGRADIENT FOR PARAMETER SET TO ZERO 16,E20.8) 
— | co,10 10 eee 
11 IF(RF.LT. UPPER) GO TO 10 
— DP(1)=0.40#DP(1) a 
ITERM( I) =ITERM(I)+1 
“— TFCITTERM(T)2LT.9) GO TO 12 
WRITE (6,2002) I,DP(1) 
2002 FORMAT (35HOGRADIENT FOR PARAMETER SET TO 10.2 I6, £€20.8) | 
FP(I)=F+e 10. O#DP(I)*PAR(T} 
mee .O CONTINUE ~~. ~~ 0 Se 


WRITE(6, 2000) (ITERLL); ITERM( 1), I=1,NP) 
O00 FORMAT (30HOBERIVATIVE MOVES##PLUS, MINUS/( LOX,8(T1,1H»IL,11X))) °° © 
WRITE(6, 700) (OPC I)» T=k,yNP) 
el =) st eee eee 
4 G(T )=(FP(L)-FI/(OP(T)#PAR(T)) 
— CALL ICLOCK(LPTIME) 
| IF(LPTIME.GT.ITIME) GO TO 50 
T JF{MFLAG-3) 60,51,60 a 
50 WRITE(6,1000) 
1000 FORMAT(28HO###se2 TIME EXCEEDED $$$$$$ | 0 
WRITE(6,1001) (P(2,1), IT=lh~eNP) 
~ LOOL FORMAT(32HORESULTS OF FIT UP TO THIS POINT/3HOX=LP8E14.5/ 
1 (3HO 8E14.5)) 
Me )0SCSC«S«éMRTTTE(G,1003P-.~.~— 
1003 FORMAT(13HOERROR MATRIX) 
- Deo I=1,NP | : 
6 WRITE(G,1OO4) (H( I,J), J=1,NP) 
1004 FORMAT(LHOLPBE1L4.5/{ LHOBE1L4.5) ) 
PUNCH 1002,(P(2,5I1), T=1l,NP) 
M1002 FORMAT(6E12.5). ae 
oo 5 J=1,NP 
_ PPNCH 1002+ (HOdel dws T=J,NP). 
® 51 PUNCH 200,(OP(1), I=1,NP) 
Mm CALL THEQRY(1,2,2,0) 
WRITE( 6,900) 
~ ANE=NE 
00 30 I=1l,NE 
m0 DEV(I)=S(I)-EC(T) 
STAN=0.0 
mepO 31 f=1,NE 
31 STAN=STAN¢(DEV(I))##2 
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on ee nm | ~~ - ~~ we oe -- ere a = - - - A ee A Re www eee em te ~~ -~ - __—- ~~ 


BA 


STAN=SORT(STAN/ANE) 
WRITE(G,9OLIISCI) ECT) OEVI(I) RIE) VCT), L=1,NE) 
901 FORMAT (2F20.3,FllelyF9el sy 15XAG6) 
~ WRITE(6,700) (O0P( 1), L=1,NP) 

ANODEG=NUDEG 

CHI2=2.0#F/ANODEG 

WRITF(6,903) NODEG ,CHI2,STAN 
“903 FORMAT(30HONUMBER OF DEGREES OF FREEDOM= 16/ 
- L34HOCHI SQUARED / DEGREES OF FREEDOM= 1PE14.7/ 
~"235HORMS DEVIATION FOR EACH LEVEL= +OR- OPF6.2) 0° 7 7 





DO 55 J=l,NP 
. DO 55 J=I,NP 
AH([,.-JI=H(1,J) 


-_—- - 


MS AH(J,1)=AH(I,J) See ee 
CALL MATINV(AH,NP,BH,0,DELTA) 
: Sewer te(6,1005) DEWTA  } — a 
1005 FORMAT({ 7HODELTA= E20.8) 
Mm  JF(LPTIME.GT.ITIME) CALL EXIT 7 7 77 OO OE 
60 RETURN 
_ te eB i ee 
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MaeORY 4 .«.-. 
MalS SUBROUTINE CALCULATES THE ELECTRONIC ENERGY LEVELS OF ARS 
EARTH TON IN NTERM TERMS OF THE OPTICAL SPECTRA. THE MAX J=15/2. 
~ HIGHER J MAY BE HANDLED IF THE APPROPRIATE DIMENSION STATEMENTS 
ARE CHANGED. THE LEVELS ARE CALCULATED FROM GIVEN CEF PARAMETERS. 
NTERM MAX = 20, MAX NO. OF LEVELS = 100, ANY POINT SYMMETRY THAT 
USES ONLY CEF PARAMETERS WITH EVEN N AND BM. 
— INPUT DATA “a 
NTERM=NO~. OF OPTICAL LSJ TERMS 
NJ=NO. OF DIFFERENT J VALUES IN TERMS 0” 
XJ=MIN. J VALUE - 1.0 
IW=DUMMY VARIABLE . in 
IOEG=1 FOR INTEGRAL Jy 2 FOR HALF INTEGRAL J 
DELTA=LIMIYT FOR OFF-DIAGONAL ELEMENTS IN EIGENVALUE SUB- = "so 758 
ROUTINF EIGENH 
NROM=STARTING POINT FOR RANDOM NUMBER FUNCTION ROM 
SVrerRY=C2 OReC3H AS The CASE MAY BE 
mots YM=C3H [ree riiiselsmiie eCaGe we =e = 5 : 
THETA(N,MI=REDUCED MATRIX ELEMENTS, N I[S ORDERZM IS TERM 


—-— ee a i ee 


— oe ee ce oe ee we eee 


— IDENTIFICATION i ae 
e AJ(1)=J VALUE FOR THE ITH TERM | | 
ee 0—~*~*C*CNNOWSVRR(M,N)=l TE THE NTH LEVEL 1S NOT OBSERVED AND ZERG 00 


OTHERWISE 
— NORMAL(M)=NO. OF THE LEVEL IN THE MTH TERM WHICH IS TO BE 
NORMALIZED TO ZERO 
— NOP=NUO. OF CARDS WITH GPERATOR EQUIVALENT MATRIX ELEMENTS FOR ~ 
EACH J VALUE DIVIDED BY 3 
HEAD(K,1)=LABELING FOR EIGENFUNCTIONS ~ 
FMT(K,y,J)=FORMAT FOR EIGENFUNCTIONS 
OA, OB, UC(K,yI,J)=1JS OPERATOR EQUIVALENT MATRIX ELEMENE FOR 
KTH J VALUE FOR ALPHA, BETA, GAMMA RESPECTIVELY, I AND J 
GO FROM 1 TO 2JN+1l, 1 IS INDEX FOR MJ=-JY AND 2JN+1l IS INDEX 
FOR MJ=J 
Mere ROUTINE THEORY {19,J0,KQ,10P) ~ ~~ | _ 
DIMENSION OA(8,16,16),08(8,16,16),0C(8,16,16),THETA(3,20) »W(16) 
DIMENSION CEF( 1641652) »GIVE( 300) »ANM( 394-2) 7AN(20),S5SIGN(2) 
BUMENSION E(100),R(100),V(100) ,P(3,40),X(20) ,QS0(3),DP(20),St tage 
DIMENSION DEF(16%,16,2) »W1(16) 
DIMENSION HEAD(8,33) eFMT( 8,12) sFMTA(12),EVECI(16) ,EVECRI16) 
~ DIMENSION C(l6,16) °}»}» ~~ i ae (ia 
DIMENSION NOBSVR(20,20),sNORMAL(20),$51(120) 
COMMON/COMTHE/NFORM,P ee Ulm 
y COMMON/COMQSQ/NE,E+S oR 
_— COMMON/COMWAV/C1 (256) 
COMMON DUMMY(281), NP 
" COMPLEX C,CTEMP,Cl. 
COMPLEX CSORT,CONJG 
CALL OVERFL(J123) 
10 LF(NFORM) 1.142 
1 READ (5,100)NTERM»NJyXJe 1 Wy IDEG,DELTAgNROM 
100 FORMAT (212,F10.5,212/E20.5513) 
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NROM=NROM¢ 1 

DO 91 [=1,NRDM 
91 IRDM=RDM(XRDM) 
Mm =6©3READ(5,73) SYMTRY,C3HSYM | 
73 FORMAT(12A6) 

DO 999 M=1,NTERM 

Mend (5,101) THETAIN,M) »N=A1l-3),AJ(M} 
mero) FORMAT(3E20.5,F10.5) 

READ(5,99) (NOBSVR(MsN), N=1516)sNORMAL(M) 


—- 


™ 99 FORMAT(1611,12)° - 


999 CONTINUE 


—— — — ee ee —_— om we ee ee —-— —— im a oe ee omy wo 


DO 1LOOO I=1,8 
DC 1000 J=1,16 
00 1000 K=1,16 
OAC IT ,»J,K)=0. 
meee OH (I,J ,K}=0. 

wetl,J,K)=0. 
(“1000 CONTINUE 

DO 998 K=l1,NJ 


—_—— _-—— < 


~~~ READ(5,103) NOP, (HEAD(Ky I) ¢1=1,33)9( FMT (Ky J) yJ=1412) 


103 FORMAT(12/20A4/13A4/12A6) 
eno 997 L=1,NOP © 


— i we we ee 


ame oe aaa ee Ne 


emia 102) 1,5 COA UKM by ly Jel O8(K oI) lig dl OC( Kol ee 


Mm 102 FORMAT (212, E20.8). 
997 CONTINUE 
_— n= K SE Nee 
LMAX=AK+XJ+0.52 
me 00 996 L=lysLMAX 
MIN=L41 
— JMAX=2.08#(AK+XJ)+1.01 
MAX=JMAX41-—L 
~ PO 996 N=MIN,MAX 
NA=JMAX+41—-N 
RermKyMAXsNI=OAIK,NA, LE) 
OB(K »MAX,N)=OB(K,NAgL) 
OC(K »~MAX,N)=OC(K,NA,L) 
996 CONTINUE 
: Heeo95 L=l,JMAX 
7 Boee50 N=1,L 
MK ,N,yL)=OACK,LN) 
Wek =N,L)=OB8(K,LeN) 
OC(K,N,LI=OC(K,L,N) 
850 CONTINUE 
995 CONTINUE 
998 CONTINUE 
_— SETTING UP OF ANM 
mo LOO). 1=1,3 
wee TOOT d=l,4 
bor LOoO)] K=),e 
ANM(I,J,K)=0.- 
1001 CONTINUE _ 
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ey 1 A=0 
i= | ee ee | 
iSA=1 
——1cikR=-0 
DO 994 N=1,3 
_— MAX=N¢1l S22 7 > — 


00 994 M=l,MAX 
meg 94 L=1,2  * 
IF(L-N) LLl,»,i1ll e994 


ee eee wee oe eee 












IF (M-L)994,993, ae 
993 ICTR=ICTRl © 
PeCSYMTIRY.~EQ.C3HSYM) 


5 [A=IA+1l 
~~ TFC TA-1Q)992,991,992 
992 ANM(N,MyL)=P(KQs 1A) 
mGea10 994  — 
991 ANM(N»MsL)=P (JQ, TA) 
~ 994 CONTINUE © 


IF (NFORM) 112,113,112 


3000 FORMAT (SHOANML/(6E20.8)) 
mer2 CONTINUE © 
me) CUS EGN(1)=1.0 
SIGN(2)=-1.0 
~ 00 900 T=1,NTERM 
PJ=AJ(1)-XJ+ 01 


IF (SYMTRY. NE. C3HSYM) 


_— oe me ee: eee ee = 


DO 903 J=1,JMAX 


113 WRITE(6,3000)(((ANM(IN,M,L), L=le2), M=le4), N= 


1 GO TO (54994, 549944994 49944994454 9944994 994 4994, 5,994)» 
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CALCULATE CEF MATRIX ELEMENTS 


Se nine <a om mee A Ce em ma et me) Sem me Sey ee ee ce ee ey i ee ee ee ce et te Oe ee ee ee ee 


1A= NP-NTERM 
Meee Gl. 1TA+1) GO TO Ge vr 


~ DO 903 K=1l_,JMAX 


{ 00 903 L=1,2 

= CEF(J,K,L)=0.0 

903 CUNTINUE 

Mme 0 901 J=leJMAX 
DO 901 K=1,J 

TF (J-K)800,801,800 





GO TO 810 
800 IF (J-K-2)901,802,803 
802 M=2 fa ee. 


SO TO 810 
— 803 IF(J-K-4)901,804,805 © 
804 M=3 
- GO -T0°810 a ase 
805 IF(J-K-6)901,806,901 
806 M=4 ie 


810 CONTINUE 
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DO 899 L=1,2 

CEF(CUsKyLI=SIGN(L)I #(THETAC I, T) *OACT I,J eK) #ANM(1LSM,LI*+THETA( 2,1) 

bait OBC IT S»>JeK) #ANM(2aM,ye LI+THETA(321)*#OC( IJ, yd KI)SANM(3, Me Lae 

Beerik, J»zL)=SIGN(O)SCEF(J-K,L) ~ _ 


99 CONTINUE 
O1 CONTINUE 


DIAGONALIZE THE HERMITIAN: MATRIX CEF 
TFC TOEG-2)500,501,500 


501 JMAX= JMAX/IDEG 


—ono-)| A ————— 
DO 503 J=lyJMAX 
KAB=1 os a ee 


DO 504 K=1,J 

moomo05 L=l1,2 co Bete ee 
DEF(J,Ke~L)=CEF(JAB,KAB,L ) 

tet, J p-LY=SIGNCL)*#DEF(CI,K yl p77 eee 


505 CONTINUE 


. —— — 


Meee AB+ O00 EEE EE EE EEL IEE EEE EL EEE ESS 


204 ONT INUE 


Ane 2 nn een ene = ++ -- 


503 CONTINUE 


| 


. 


me te re 
ee 


49 


50 


es 


-_ ~ 


Seen en 


50 
50 


MeemroraAXeEQ.1) WCL)=DEF(1,1,1) 22 ° DOO 

IFC JMAX.EQ.1) GO TO 506 
Meee cIMAX.NE.2) GO TO 499 ee _ 

TEMP= SORT((DEF(1, lL» 1)-DEF (2, 2ylilus2er4. OSIOEF U2, l» gl lees EF (2, yee 

—Laee2)) 

W(1)=(DEF(L,1,1)+0EF(2,2,1)+TEMP)/2.0 

W(2)=(0FF(1,1-1)4+0EF(2,2,1)-TEMP)/2.0 — 

GO TO 506 
9 CALL EIGENH(DEFyWeJMAX,DELTA) EE eee 

coe1O 506 
O MZ=0 icy ie as eae em Toa * a 

DO 515 M=1,2 

Meee = ( IJMAX—-1L)S 2000 ae 

IF (M.FQ.2) JMAX1L=(UMAX4#1)/2 
Be oo ae , eee ae ee 

Mie EQ. 2) JZ=1 

00 507 J=1lsJMAX1 © . i i 

KZ=? 

Memeo 2) KL = p00 a ee 

DO 508 K=l,J 

DO 509 L=1,2 

DEF(J,K,LI=CEF(UZ,KZ,L) 
9 DEF(K,JeLI=HSTIGN(L) @D0EF (J eKeL) 
8 KZ=K2+2 


Mm 2= S242 000 i — - 


meet ec)) 


[FC JMAX1.2EQ.1L) WIC lL)=OEF(1,1,1) 

Pet JIMAX1.EQ.1) GO TO 511 

IF (CJMAX1L.NE~2) GO TO 512 
TEMP=SQRT((OEFCLe Ly LI-OEF( 222,11) )8#2+4.08(DEF(2e1leli) ##2*DEFl2y1le2) 


oe er ee ee | -_—— E —— ee . - —_—_—— 





ewww ee ee _ — —— oe ee te oe - = ——_ _ RR tn - -— — —  — — 


B= 


etl) =(DEF(1,1,1)+DEF(2,2,1)+TEMP)/2.0 
W1(2)=(DEF (CL ,1,1)*DEF(2,2,1)-TEMP)/2.0 












Sea 10) 5Sl1l 
Beers CALL EIGENH(DEF,WlsJMAXI, DELTA) £  ......©=»©6©m6—= a 7 
Srl DO S10 J=zl,JMAX1 
_ MZ7=MZ+4] ee . oo a 
510 W(MZ)=W1 (J) 
™5)5 CONTINUE a re + a — 
S06 KMAX=JMAX-l 
Mee, K= pA 
KPL=K+l 
eng 17 J=KP],JMAX = = °°”. OT a 


IF (WOK )-W(J))17,17,18 

TTS =W (Ke ee eS eee 
W(K)=W(J) 

evigj= 'EMP 

Mm i? CONTINUE 

— TNORM=NORMAL(L) 

TEMP=W(INORM) 

D0 69 K=l,JMAX 

69 W(K)=W(K)~TEMP 


ee ee ee ee ee ee es a ee me me i ee me ae es ee ee ee ee ee ee es ee ea ee ne ee a 


en eee ee cee ei Se ee ee ee ee et ee ee em re ee ee rt ee re ee ee ee i ee ee ee ee et ee wee es ee es --/- 


™ OO 50 K=1,dMAX a = 

: IF(NOBSVR(I,K).EQ-1) GO TO 50 

Mee) 6 6StiS)=WwiK) ee Perak Baie a Fe SS ae 
mt S)=SC(1S) 


IS=1S+1 

50 CONTINUE 

eco 10° 53 

me 52 IF(IDEG.EQ.2) JMAX=JMAX/2 

Mm 53 1ZS=KQ a ee ee 
IZT=IA+I 

Meter Q.EQ.IZT) I1ZS=J0 
DO 51 K=1,JMAX 

he LTFt(noBSvR(1,K).EQ.1) GO TO SI 
SrSA)=SL( ISA) *#PL(IZS,I1ZI)} 
PsA=1SA+) : 

51 CONTINUE 

mo )60CSC«ir FC TQOP.NE.O) GO TO 900 


A Se a S-Ni 
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j [F(IOEG.EQ.2) JMAX=2*#IMAX 
, ~ CALL EIGENH(CEF»>WeJMAX DELTA) 
00 6005 J=lyJMAX 
> 00 6005 K=1,JMAX_ 0 ee - 
[INDE X=J+K #JMAX-—JMAX 
mO6005 C(K,J)=CLUINDEX) — se 7 : : 
6006 KMAX=JMAX=-1 
C€  &MULTIPLY C BY PROPER PHASE FACTOR ee 
her 46000 J=1,JMAX 
, TEMP=REAL(C(1,J)) cw a 
IF(TEMP.EQ.0.0) GO TO 4001 
~— CTEMP=CSQRT(CONJG(C ( UMAX2J))/EC(1,J)) 
GG, 7O 4002 a a | ae 





















——— oi 


KPL=K+ 1 


— ow ee 


RETURN 
~~ END 


ih ee wee we ee ee eee ee ee ee ee eee ee es ee ee en ee eee eee ee mew eee eee i + —— 


a 
eee 
—— 


900 CONTINUE 


-— 


00 6017 J=KP1l,JMAX 
[F(W(K).-LE-W(J)) GO TO 6017 
TEMP=W(K) 
~ W(K)=HW(S) 
W(J)=TEMP 
mee o018 L=)]"JIMAX 
CTEMP=C(L,K) 
mms KK) =C(L, J) 
6018 CllL,J)=CTEMP 
6017 CONTINUE ~~ 
WRITE(6,600Q) (HEAD(IJeK) »-K=1,33) 
6000 FORMAT(33A4) 
Homo00] K=1,12 
mOOl FMTA(K)=FMT(1[J,K)~ . 
HoeoOO2 J=1,.JMAX 
~~ DO 6003 K=1l,JMAX 
EVECR(K)=REAL(C(K,J)) 
6003 EVECI(K)=AIMAG(C(K,J)) 





—_— ee me 


~ —  e —- ee e  e 


ee 


ee 


—_ 


ace 


4001 CTEMP=CSQRT(-CONJG(C(KMAX,J))/C(2,J)) 
4002 DO 4003 K=l,JMAX . 
4003 C(K,J)=CTEMP#C(KyJ) 

4000 CONTINUE i — 
DO 6017 K=lyKMAX 


ee ee 
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eee ee 
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_ — ~— ~ ~ _ _—_ — 
( 
mm meee ee ee ee ee ee ee er ee ee ee eee — -_—- 
-—- — ~ teem eee eee eee ee ~— ee ee ws oe — mee ee ee ae ee ee me meee -_——- —— — 
ee ee meee ee ee ee _—— —_ — - -- ——_— = —_ —-<«-. - — rr -—s oo 


le a er re rr ooo ten — — ——_— — 
_ — _ 
— . een a ee oe 
ee en ee er ———e — 
ee ~— I —_— 
— ee — — a —— rr _ 
~ ger pep — ~ee 
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dit mm ee em i et eR eS ee ee i ew a eee ee re ee ee ~~ ——— ae ~—<—- ——— -_ -— ar wee 


| BSCUAR..«. 

— FUNCTION TO CALCULATE VALUE OF Q#e#2, WHERE CHI ##2 [S MINIMUM VALUE 
OF THIS FUNCTION 

~ FUNCTION QSQUAR(1Q,JQ,KQ) 
DIMENSION E(100),R(100),V( 100) »P (3,20) X(20) ,QSQ(3)-¢0P(20),S8(1500) 
DIMENSION GIVE(300) | . a 7 “a te” 
COMMON GIVE,S,EeReVeP eo Xe QSQsNE eNP91Q,JQ,KQ,10P,XTEST 

~QSQUAR = 0.0 vs, hg ie 
CALL THEORY(1Q,JQ;KQ,1) 


eT. l=1,\c¢ — °°» °°» © 


C 


ee ae ee ote 





a = a 
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a te ey ee ee a re ee En Se re SR te ee a SS SE dy RR RS oe Sa ci te Se rds eee ey ee ee Ss i eg ee ee ee ee ee ee ets em ey ee —— meee ee 
ee ee ee ee et kk, ke te dd ee ee eh at ms = _—— 
[<a are ee Ger es ae es ee ee es es ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee re ee oe ee ee ee ewe is ee ee ee ee re es ec ee eee ee ee ee a ee ee ee ee eee ee _———- 
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——— _—_—_ _ ee eee —— eo —- a aia ee en cd eT ee A ars ce esas ne ee —- ee ee roe — oe ~~ — a tees -_ - _ 
SO eee te ee ee eee ee ee a ee ee — ee ee ee ee ee ee _—— oe ome _ ey —_——— eee _ ——_— — - _ _ ~ 
Oe ae ey et es rm en ee eee ee eee wee ee ee ee _—s ee wee ee ei ie _ me ee wee ee -_o - - ~ —_— - oe -_ a -——_— eo on + 
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ae . on ee 
BEIGEN... 
THIS IS A LINK BETWEEN A PROGRAM CALLING EIGENH AND HERM. 
SUBROUTINE EIGENH(CEF,W,JMAX,DELTA) — e a 


~~ COMMON/COMHER/B(16),H( 200) 
COMMON/COMWAV/C 
COMPLEX H,C,CMPLX 
DIMENSION €(16,16),CEF(16,1622),W( 16) 
_ IREGIN=16-JMAX eae 
Ooo) i l=l1,JMAX 
Meme 1+! BEGIN EE ee 
Mert S )=CEF(I,1,1) 
hey : 
IMAX=JMAX-1 
mm  '00 2 J=1,I1MAX a as reer : % 
JPL=J+1 
eee K=JP1,JMAX 
fer =CMPLXICEF(J,K,1),CEF(JsK, 2)) 
ett LN ee ae = ee oe 
KZE=I+1/2-1 


- — ten — - - et 


—_—_— - 


a emt areca mK | APU ag em mm 


_— Meee =(0.0,0.0) “a © Tee 7" oa 
pmeeZE+1)=(0.0,0.0) 
~ CALL HERM(H, JMAX,C,0,D0ELTA,zIT) _ 
mest =1, IMAX 
MITT OC CTN ee 
gett) =8(15S) 
me RETURN  # © OCCCCTTTT a Ss =e . a 
END 


_— — a eae — eo 2 EE oy EEE Ee SS ee ee 





m3 







HERM.eee 

* THIS SUBROUTINE DIAGONALIZES A HERMITIAN MATRIX AND OBTAINS ALL 

: EIGENVALUES AND EIGENVECTORS. IT IS SHARE NO. 884 PK HMEE WHICH 

~ HAS BEEN MODIFIED SO THAT IT MAY BE CALLED FROM A FORTRAN PROGRAM | 
ONLY THE MODIFICATIONS ARE SHOWN HERE. 

~ THE CALL STATEMENT TO BE USED IS CALL HERM(HsN,U,PR,DELTA,IT) 
ALL ARGUMENTS HAVE SAME MEANING AS DESCRIBED IN SHARE NO. 884 

~ WRITE UP, PR SHOULD ALWAYS BE ZERO AND IT IS THE NUMBER OF ~ 

ITERATIONS 

MS AVE 1g Do ee 


_—_— 


SXA XR4;54 
—_— eo, en Fee ee 
STA Hl 
— momen G4, 4 0 ee ee eee me an 
LGL 18 
fees t> Hi 
le » 4 
MU iD 0L0L——————— A Oe 
CAL *® 654 
—_— MMe OE ee 
STD H2 
Memenim 7,4 
SLW H3 
MSO HIMEE»40 00 ae 
Hl PZE #*,,88 
-m  +H2 UG ES 7 “oe 
H3 &SS 1 
TERT 6 
CAL COMMON+20 | 
fg Beesl We B,4 Pope © ee “ - a 
RETURN HERM 
— FIRST CARD OF SHARE NO. 884 FOLLOWS THIS CARD  — 7 = 
* LAST CARD OF SHARE NO. 884 PRECEDES THIS CARD 
COMMON BSS 23  }  — —_ a ee ere : ; 
END 


_-— _ = te ee ee ee ee ee ee ee er re oe we eee ee a —= ee eee me, me Se Se tee ew ee ee es ee ee ee ee _—— -_—_—— ~~ 





The required card decks are: 
Gee 
RME 
SIX] 
Died: 


PAC 


The programs in decks 3 and 4 were written by B.A. Zimmerman. 


The listings of decks 1-5 follow. 








. —~=— = 2 oe 








an ee 


_ awe oe 


1 


ma a 


= — ee = ao 


E 


= 


em ee ee 


ine 


ee INPUT DATA 


See CU «C2 4, 


— ee ee 
, ee ie Ge 
100 - 
ee eee 
wwe we ow ee oe oe 


1000 — 


WRITE 


te ee ee ee ee ee 


-95- 


ICA 
THIS ITS A PROGRAM TO CALCULATE THE OPERATOR EQUIVALENT FACTORS 

IN THE INTERMEDIATE COUPLING APPROXIMATION 

FOR TM IV OR PR IV GIVEN THE SLATER INTEGALSs SPIN-ORBIT = 7 ~~ 
COUPLING PARAMETER, AND PURE LU-S OF FACTORS. 

REF. FeH. SPEDDING, PHY.REV. 58,255(1940} AND GRUBER AND CONWAY | 
Je CHEM. PHY. 32, 1178(1960). 


Ely £2, E3=RACAH'S LINEAR COMBINATION OF SLATER INTEGRALS 
‘6 a 

ZETAL=SPIN-ORBIT COUPLING PARAMETER 
DIMENSION AS(3),AL(3) i 
DIMENSION S(343)5 EV(393), E(3) | 
DIMENSION ERASS (20)4A(3),8(3),G(3),AE(3),BE(3),GE(3) 
READ(5,100) E1,E2,E3,ZETAI 
Barats F2009) ae 
F2=(E1+143.0*E£2+11.0#E3)/42.0 
F4=(6.0*F2-39.0*E2-E3)/11.0 
F6=(7.08F2-18.0"E2-3.08£3)/77.0 
MeyaA=-ZETAT/2<0°-°- 
WRITE(6,1000) El, E2,E3sZETALsF2,F4:F6yZETA 
FORMAT (4HIEL=LPE14. 7, 5X3HE2=E14.75 5X 3HES=E14. 7, SX6HZETAT=EL4. 7/7 
14HOF2=E14. 7, 5X 3HF4=E14. 7, 5X 3HFO=E 14s 7 y SXSHZETA=E14.7) 
PURE ELECTROSTATIC ENERGY LEVELS © 


a a ee ep ee - — et eee eee eee —e ee ee ee 


- A mm a me ee tee nee ee a me em ee ee ee ee ee ee ee ee ee ae ee en ee es ee ee ee ee 


_ me ae ee ee ee ee mm ee me ee ee eee ee 


—_— A a a ee -—_— = = 


E3H = 0.0 

apes 15.0 ® F2 + 18.0 * F4 — 273.0 # FG  °  &°& ae 
BlG = -— 5.0 * F2 + 148.0 # F4 + 91.0 # F6 
El1D = 44.0 # F2 — 48.0 * F4 + 728.0 * F6 
Ell = 50.0 *® F2 + 60.0 #® F4& + 1460 # F6 

meme = 70.0 * F2 + 84.0 #-F4 = 1274.0 * FO ° # ©) ae : 
ElS = 85.0 «* F2 + 249.0 # F4 + 1729.0 # F6 


~COULOMBIC AND SPIN-ORBIT MATRICES, CHARACTERIZED BY J 
Situle) ) E3F + 4.0 # ZETA 
S02 | Ea ae = 
3,3) E3P - ZETA 

oe, 1) 2-0 * 2.449490 * ZETA 
0.0 


3.0 9 eal 1s ae EA 


a 
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ee ee ee ee ee ee ee ee 


omen a naan 





ea ei a a 





nonuiuwouiw ea 
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ee ee ee ee os ee ee re me es er es + 


| | | | 


mam hw m= © — 


l 
(6,400) 


a a A ee ter gaye 


ee ee eee ee oe ee - co eee ee ee ee ee ee 


ee eee ee ee ce a ee ewe ee ean ee ene ee a ee ee ey ee ee ee ee oe oe oe 
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20 Ge 


meeOO FORMAT ( LHI 12x GHENERGY SISA 3SHSF2 15X 3HID2 15% 3H3P2 0137 


ISHALPHA 14X 4HBETA 13X SHGAMMA ) 








Co 10 30 

Memesti,l) = E3H + 680 # ZETA ©" © © © OOS 
Siz2,2) = £16 

“"§(3,3) = E3F - 3.0 # ZETA : 
mms.) = 2.0. .* 3.162278 / 1.732051 « ZETA 

me S(3,1) = 0.0 <2 

$(3,2) = ~— 2.0 # 3.316625 / 1.732051 # ZETA 

mum. 2) = S(2,l) °° °° a 
Bl» 3) = Os 0 
Bites) = Sta “Sa - ar : 
pe 1} =5.0 
met2)=4.0 r. ie eng 
AL(3)=3.0 

4.0 # | a SS 
NORDER = 3 

UPL AG = 2 ee err a ae 


~ 500 FORMAT (1LHO 1L2X6HENERGY 15X 3H3H4 15X3H1G4 15X 3H3F4 13X 


em 


ferme crepe pm tm — - yee 





——_ 2 a oe 


ee ee 


WRITE (6,500) 


LSHALPHA 14X 4HBETA 13X SHGAMMA  ) 
ec 1o 30° _ 
PO Stl, 1) 


ome, 2) 


Ell 
E3H -~ 5.0 # ZETA 
- 2.449490 « ZETA 

P Siow. & - 


ee ee ee eee ee ee eee ee re ee ee ewe ee es — ee ee _- —_— - — 
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behead 


© 
e 
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[FLAG = 3 
WRITE (6,600) 


600 FORMAT (1LHO 12X 6HENERGY 15X 3H116 15X 3H3H6 13x 


LTSHALPHA 14X GHBETA 13X SHGAMMA ) 
TEMP=SQRT((S(151)4S(292) 8 #24+4.08(S5(152) 0 #2-4.08S (lel )eS(2,2)) 
mem) =(5(1,1)4+5(2,2)47EMP)/2.0 
Be =(S(1,1)¢S(2,2)-TEMP 172.0 


-_ 


Meee {=1,2 
MeMe=(E(1)-S(1,1))/7S5¢1.2) 
EV(L,T)=1L-O/SORT(1-0+TEMP#82) 

12 EV(2,1)=TEMP/SQRT(1L.O+TEMP##2) 
mo TO 31 

__30 CALL EIGVV(S,EV,E,NORDER,ERASS) 

Sm 31 DO 25 I=1l,NORDER “- aT — 
BUT V=RMECAS(I),ALCIYFALUID, AJ, Ads 2.0) 
POLI =RMECASTI) pAL(1 }sALT 1) 4Adsh),4-0) 

memGt(l)=RME(AS(I),AL(1),AL(1) -AJ,AJ,6-0) 
IF (NORDER.EQ.2) GU TO 5 

_ ALBSRME(AS( 1) yAL (1) AL( 3), Ade AS10260), 


a eee me ee ee ee ee eee - ~ é —— ee ee ee em ae 


—— a we ew te ee ww ee 





Gy 


5 CONTINUE 





300 FORMAT(E20.8) 


WRITE( 6, 300) (ACT) BCT) 9G (LNs ~~ T=l,NORDER),A1L3+613,613 





BL3=RME(AS(1),AL(01),AL(3),AJ»,AJ,4~-0) 
GL3=RME(AS(1)sAL(61)2AL(3),AJ,AJ26.0) © 





— NS pp ar ma 


_— DO 3 I = 1» NORDER  }}  — 
metl) = 0.0 
me)6CMUBME( I) 6 = 0.0 -_——s 
fet) = 0.0 
Seoes J = 1, NORDER 7 
Mer) = AE(T) + AlJ) # L[EV(J,1))*82 
~~ BE(I) = BECI) + BIS) #® (EV(U,1))##2 
PeGeii) = 


we ee oe 


PeeNORDER.EQ.2) GO TO 3 


AeE({) = AE{I) + A13 © EVEL, 1) 

“ BE({) = BEC!) + 613 #® EV(1,!1) 

etl) = GE(I1) + G13 * EVI1,1) 
-_ 3 CONTINUE — oe 





IF (NORDER.EQ.2) GO TO 63 
WRITE (6,200) (E( I), (EVCJ,I) 44 
1 { = 1LsNORDER) 
200 FORMAT (1LHOLP7E18.7) 
eomro (10,20), .IFLAG 
63 WRITE (6,201) (CECI), (EViU,T Ye 
1 I = 1,NORDER) 
201 FORMAT (1HO1P6E18.7) 
GO TO l 


| cite 
END 


in ee ain es 


a) ——_— — 2 - 


— ee eee ee ee eee ~—— — — 


ee ee ee we ee es a — ee = 


GE(T) + GlJ) # (EVI0,1)) #82 


— a ee ee oe ~— — ene eae — oe we _ ee eee oe ee — ~-_~ 


# EV(3,1)2#2..0 
*® €V(3,7)#2.0 
* EV(311)e2.0 


—m ee -— oe ie ~ — one ——_ —-_- — - ~— 


L»NORDER), AE(I), BECI), GE(T), | 


a 


LyNOROER), AE(T), BECT), GEtI), © 


me ees ee — = ee ee ee we es ee ee we ee 


en a ee eee ee ee ee ee ee es se re ee eee _——e— oe A oe oe fe ew Se ee ee ee ee ee ee ee ei ee ee ee eee ee ee ee ee ~ -_-_ = 
| 
; a oe ee ew me we we ee es eg ee em ee ee ee re ee er es we we ee ee ee ee ee ee es ee ee ees es ee + Se ee we ee we ee ew ee ee -_— — _—— — —_ 
ee: Se ee 
jw ewer ww ww me wee we we ee eee we ee — oom - -_ m= - _ -_— = — _ — —_—— = =_ —oo— 
a tte eee ee -_—_—— —=— —_ = — 
a ee ee a er te 
me cere ce ee we eee ee ee ee ee we we ee we a ee ee ee ee eee — ee ee; Cit} Se ede te se es eee See ooo we os a oe oo — —_ ee ime —_ 

-—— oo mem ee me me ee re ee ee ee ee ee ee ee _ —_— = -_— ~ - me oe ee ose _— —_ 
a a ae 
_ ~ _ ~_ ~ ~ _ oe ~_ ~_—— os eee ee Re ~-- — — _ 

— — = ~ — -_——— _—— — _ - 

eee 0 a i er ee mm — 
= — we ee a _ — — — oe _ —~ ——_— ~ 
—— = _~ ee re a ee EO Ee ee ee PE EE ee Se ee ee EE ee ee ee ee ee ee ee ee ee ee ee -— = = _— 
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Pas Ee ee 

THIS IS A FUNCTION TO CALCULATE THE REOUCED MATRIX ELEMENTS UR 
OPERATOR EQUIVALENT FACTORS FOR THE 4Fl2 ELECTRON CONFIGURATION. 
Bedi ELEMENTS ARE OF THE FORM CS ey J LINE Do Spee yp ys ~ 
REF. BeRe-JUDD, PROC~ ROY. SOC~ A2415414(1957). 

BeeeraisS VERSION [S FOR J=J*. ‘ 

FUNCTION RMECAS,AL,BL+AJ,BJ, AN) 

meeSC TEMP U=SQRT((02.08AN+1.0)4(2.08#B5t1.0) ) 

TEMPL=SQRT((2.0#AL+41. O)e(2. -O*BL tl. -0)) 


ee l tn et er ee ee ete ee ON i tn mre ne em a mere rt ep ee ee ee 





Pe= SIX JCAL,AJ,;,BL»~BJ,ASsAN) #(—-120) © (LIF+L24+dle52) 
Meee SIX IC 320, AL ¢ 3200, Bly 300, ANI OU-1 Ole (Cl+lLetsJiei2) "Se 
N=AN 
_ J=AJ pat : 7 —— ees . ~~) ae 
Bere T=SORT( FACT (24#J-N) /FACT( 28 J4+N+1)) 
RMEL=(-1L.O)#e(IS—- J+t2 )#TEMPJ#®TEMPL#SFACT#WwleW2 
GO TO (lr2ele4e156)5N 
2 RME=16.0#SQRT(7.0/15.60)#RMEL 
RETURN 
4 RME=-32.08SQRT(14.0/11.0)*#RMEL — ——— 
RETURN 
© 6 RME=1280. Dope ett )o acetone 
RETURN 
Meee te (6,100) 
100 FORMAT (32HORME DOES NOT EXIST FOR THIS N ) 
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netroots ee 


aoe ; : : a ee oe —~ re el a 
OM pn a ee - te ee re ee 
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oa Js. « 
A FUNCTION FOR SIX J SYMBOLS 2/25/64 . 
ROTENBERG ET. AL. PAGE 13 EQUATION (2.3) 
fens «INPUT SIXJ(JT,J2,L2,L1-J33;013) IN FLOATING POINT “—" wee em 
REQUIRES DELTA AND FACTORIAL ROUTINES 
FUNCTION SIXJ(A,8B,C,D;,E-F) a . "2 ; - _— 
TRI1 = At+B-E 
mmr RiL)2;st.1° | nt 
2 SIxXJ = 0.0 
Vuk  — 
1 TRILL = A-BtE 
Mee weiF(TRI1)2,3-3 °° Set ae ep er Co ee an — 
3 TRI1 = ~A+BrE 
— PE(TRI1)29454 9 9 er eee 
FIRST TRIANGULAR TEST COMPLETED 
om—A 2? =-D+C-E “ © 
Met TR12)2,5,5 
Murs? = O-C+E 7 
IF(TRI2)2,6,6 
TT OC 
BEAGTRI2)25 797 
SECOND TRIANGULAR YEST COMPLETED 
7 TRI3 = A+C+F 
me )|6)©6™SCU FC TRI3)2,8-8 
8 TRI3=A-C+F 
Mr )60SCU TIF (TRI3)2,9,9 
9 TRI3 = -A+tC+F 
Meer (TR13)2,10,10 
THIRD TRIANGULAR TEST COMPLETED 
™ 10 TRI4 = O+B-F ee; eC ee 
MeCTRI4) 2, b2y V1 
ll TRI4 = O-B8+F ers - = 
IF(TRI4)2,12,12 
Peer 4=-0+8+F °° ~~  — 
. DEC TRI4)2,13+13 
—_— FOURTH TRIANGULAR TEST COMPLETED ~~ Te : 
13. DELI DELTA( A, BoE) 
Seetyet 2 = OELTA(D,C,E) ~~ 9 ~~ es 
DEL3 DELTA(0,ByF) 
MER 4 = DELTA(A,C,F)Y) .|@©§©9”™”™©9@W©09W. ©) aa 
DELX DELL*DEL2*DEL3*#DEL4 


a oo ee ae oe ee 


ee oem ome Se Se oe 





er ee ee cree ee rm ee re ee ee ee ee ee me ee ee ee ee ee re ee er ee cm erm wre es ee ee et ee ee ee ee ee ee ae 


ee ee ee we ae eae ee ee ee 


a ee ee 


iow. wo ow tt 


Om em ee ee 


N = A+B+C+0 ee eee 
PHZ = (-1.0)*8N 
—_— SUM = 0.0 ee / 
| k= 0.0 
SS = eee 
a [F(S1)22,516,16 | 


i M = Si 

pole = FACT (M) 
S2 = C+D0-E-AK 
— TF ES 2)225185 18 


; ia ae ee sls - «=. 


ne RR tA TE A, 


Me mt ee ee ee ee re ee ee ee ee ee er ee re we ee ee es ee ee ee ee re -  e  e e e e e e  e — ee ee ee Oe ee aa ewe 


























- l 0 O - 


a= S2 

mee = FACT(M) ~ © 
Bo = AtC-F-AK 

J M153) 22,19,19- 
fen = S3 

ees s = FACTI(M) 

S4 = 0+B-F-AK 


ee me mm wm mt te + rm ee mm 





ewe ae - = > - ee ee _— . - = - - - 





ME t54) 22,20, 20 : : 
ees. = FACT(M) ———~™O 
— S5 = -A-DtE+FtAK 
Mei f(S5)17,21,21 Sees ee LL ee ee 
Pome = SS 


foo = FACT(M) 

S6 = -B-CtEt+FtAK 
Beets 6) 1,235 23 
23M = S6 


o— = = om 





me ee ee ee i ae em ee me ee ee ee a ee ee ee ee 


FS6 = FACT(M) © a 
N = AK 
Mee? = (-1.0)*#eN ~ © © = 
TOP = A+B+4+C+D41.0-AK 
a ie= 10P 
FTOP = FACT(M) 
‘TOP = SPHZ#FTOP . 
M = AK 

~ FAK = FACT(M) 

, DENOM = FAK#FSL#FS2#FS38FS4*#FS5#FS6 
~~ SUM = SUM + (TOP/DENOM) 

17 =AK = AK + 1.0 
me GO 10 24 
22 SIX = PHZ*DELX*SUM 
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a = = me — _ oe ee me ee eases ee ee ee ~ — - _—_— — — owe 


won ae OU mee eee ee ee aa as we ee ee eee ee ee ee eee ee ee 


~~ mm mem me me eee we oe ee ~ - - a — 
RETURN 
END 
a ee 
_— oe -—S = —_-~- -- - - = - ——_—-— ~ - —_ a) - xe o - — ; = = —— _ 
- -- - - — - 
es a ae IR ee ne eee ae ee ee 
ee ae ee ew -- = —~ - ms —_- = _ - -——-~+ & A 5 eae — — ——w ee ewe eee 
ee ee ee ee ee ee - ——— —— oe _——s ~ ~ — eo en eo ee 
eee ar 6 ee ee ee 
— _- -_ -_ - 
= ~ ~ a oe 
a ee mer ee tr aa PS eS — ee AR RY — ose 
ae -_ - ~~ me ee ee 
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DELTA... 
~ FUNCTION DELTA(A,B,C) © 

ROTENBERG ET. ALs PAGE 13 (2.4) 
FUNCTION DELTA(A,ByC) a ial 
m= A+ B- C 
i Lana So < ee et tt ee Se lt - ee 
FS1 = FACT(M) 
Se - A+€- 8 

S2 | 

= FACT(M) 


Meee Ut CCK A ; . — ae: —_ 
M = S 3 7 a) oa) ae "os ana ial ii moe ae ao men rae oy T° a. trai een) Tei Gn at es ll 


FS3 = FACT(M) 
~~~" DENOM = A+B+C+4+1.0 
| M = DENOM 

Meee = FACT(M) 
DELTA = SQRT((FS1*FS2#FS3)/FD) 


~~ ee en ee eee oe es ee ee em ce ee ee cee es ee ee ee ee ee een em ae ee ee 
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ee se ee ee we ee em ee ~~ = —_—_— _ — _— a _ ~- 
cme mee ree ree ee ee ee we eee ee ee ee ee ee ee ee a we oe _ — = — eee ewe ee ia eae — me ee owe — 4 ee en ee _— ~—-— = 
4 ee eed ee eee a wee ~ ee wee ee me ee ee ee —_ — — =- = 
— 
_— _- ~ ~ od ~ = - _ — ~ — —_—_— == a _ _—— = = 
—— — — = - —~ ~—— ~— — ee ee ~ = -_ 
ee a rm ne ee a ee a ee 
-— a - _ -_ “=~ = ar - - - - i _ 


- ~_—— wwe ee eee ee ee ee re re ee ee ee mo — awe a eee re re ee ee ee eee we eee —_——-= _ 





ee 
CALCULATES N FACTORIAL 
_ FUNCTION FACT(N) 
es a 
eee ac 1 =0.0 
STE A. LT.0.0) RETURN © 
me 7 =1.0 
 TF(A.EQ.0.0) RETURN. 
FACT=FACT#A 
meme) 1=1,100 
a | 
ee-A-B- 8} }=30—)?>= 
imi <.EQ.0.0) RETURN: 
Mec tT=FACTs«C 77 
1 CONTINUE 
~~ WRITE (6,100) 
LOO FORMAT (1LOHOFACTORIAL OVERFLOW ) 


—— ee ee ee re em rt re ee ee we ee ee ee eee ee ~~ _—— = 
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ee ee ee ee ce me ee me cr ee ee we ee em ee - 


mem we ee re ee we eee ee ee ee ee wee ee eee ee ae ei ae i ec ee ee me + 


ee 


- see ee -_ _ 
~~ -_ i oan 
eee ee —— — weer ee ee 
=a oe = == _—_— = -— = _ _~ — — — — —_— = 
ea | _ = - 
ae ec re mn ree ee a er ee 
s 
- 
a ee ere ere a = A I oe ay 





f\2 
-~ PUI = 
SoS 


QTAVE 


The required card decks are: 


SEARCH 

COREL 

OUTRUE 

DERIV 3 

JPEOT 

CEFMAT 

THEORS ss 

QSQUAR see listings for OPTIC 


HERM see SHARE no. 884 


The listings of decks 1-7 follow. 








SEARCH... 
THIS PROGRAM USES THE METHOD OF DIRECT SEARCH TO OBTAIN BEST 
VALUE OF CHI ##2 
me )™CSCNPUTT DATA | . EE 
NE=NO. OF ENERGY LEVELS MEASURED 
NP=NO. OF SHIELDING PARAMETERS 
NPP=NO. OF SHIELOING PARAMETERS WITH NON-ZERO STEP SIZE 
“ITIME= MAX. XEQ TIME ALLOWED ; 
E(I)=MEASURED ENERGY LEVELS : 
Meee )06CéaR CTY =ERROR “IN MEASURED ENERGY CEVELS ~~ ~~" (a 
| V(IL)=LEVEL IDENTIFICATION 
P(2,1T)=INITIAL VALUE OF PARAMETERS 
X(1)=STEP SIZE FOR EACH PARAMETER FOR INITIAL CYCLE 
XMIN=2e#N , WHERE N IS MAX NO. OF CYCLES ALLOWED I.E. FINAL 
STEP SIZE WILL BE X(1)/2"#N a 
~ OIMENSTION E(100),R(100),V( 100) ,P(3,20),X(20),Q8S0(3),0P(20),S(1500) ~~ 
DIMENSION AB(20) 
~~ DIMENSION GIVE( 300) 
COMMON GIVEsSyEeRe Ve Pe Xs QSQeNEoNPe IQs JQeKQe IOP sXTEST sNFORMyNMOVE GA 
COMMON /CGOMTIM/ITIME re ee — 
READ (5,499)NE,NP,NPP,ITIME 
SMR MATCIZ/13/13/16)  ......©°=»=£=£© 0 
CALL ICLOCKC INTIME) 
~ LTTIME=INTIME*ITIME 
mer OD «(5,800)(E(T),RII),VC(1), 1=1,NE} 
800 FORMAT(3F20.5) ~ 
READ (5,500) AA, BB, (AB(I1), I=1,NP) 
500 FORMAT(12A6) 
, NFORM=0 
MIR EAD (5,200)(P(2,1).X(1)- T=l,NP) °  .8§©)>©m(©D [m©—©—©—CU 
(200 FORMAT(2F20.5) 
Meee | 6©'READ)hlU(S,300)XMIN  } }©=~— ee es ee 
300 FORMAT(F20.5) 










wee ee eee 


— — =< ee / oo oo ee me ee ek ee ee ee ee ee ee ee ee ee em ee le - -_ ~ ~_—_ 


same se ep eee ee oe eee 


a mg 


~ NMOVE=0 — 

Moi 2)=QSQUAR(1,2,2) 

-_ NFORM=1 : as 7, ~ ae — 

| A=NE-NPP 

Mm  £CHI=QSQ(2)/A : “i 7 
WRITE (GuuOiONkGne Cl jal le Ned 

Mm 700 FORMAT( 1H] » (10X%sA6—14X% 9 Abs 14X%_y AG, 14K pAb L4X,AGV14X3AG))~ 1? 
WRITE CG cer Oly) 

701 FORMAT(9X,7H CHI *##*2//29H INITIAL VALUES OF PARAMETERS) Si 
WRITE (6,702) (P(2,1),1=1,NP) 

™ 702 FORMAT(6E20.8) 3 © a | aa 
WRITE (6,703) CHI,(X(1),{=1,NP) 


— 703 FORMAT(LPE20.7//24H INITIAL VALUES OF STEPS/{(0P6E20. 8) ) 
C FORMING PATTERN 


XTEST=1.0 
meet NZERO=0 ee eee 


NFORM=NFORM+41 
_ 0O 2 {[=1,NP 


ee en th Ra =e er te i SS RS Pe St =- —— 


—— oe we oe —— ——_ — eee we ee ee tee eS See UE egirEmE ete Se nM =e wae ee a Pee ee SS Se wee ee ee we wee ee ee ee ee we we 








= ee - ees = ee err creer rr a ee a ee ner re a ee ee eee wae lO 
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et x(1))50,5,50 
~  §0 CONTINUE ? ea = 
P(l,1)=P(2,1)-X(1) 
Mt 1)=OSQUAR(I,1ls2). 
Pa, 1)=P(2,1)+xX(1} 
fe OS50(3)=QSQUAR(I,3,2) , 
m0SQ0(2)-08SQ(1))3.3+4 
Mees 1F(QSQ(2)-QS8Q(3))5,5,6 _ - oo 
4 OP([)=-1.0 
Memeeor7)=OSOQ(LY = 
eo TO 7 
memeiP()] )=0.0 
NZERO=NZERO+1 
me GO TO 7 ee ee Les a | an 
por ti)=+1.0 
eo (2)=0S0(3) 
Meee 1)=P(2,1)+0P(1)#xC 1) 
2 CONTINUE 
CHECK TO SEE IF PATTERN IS NON-ZERO 
MANZERO—NP) 6,958 (00 eee 
9 XTEST=2.0#XTEST 
~ CALL [CLOCK(LPTIME) 
BeelPTIME.GT.~IT{ME) GO 10 i) 
oom F (XTEST—-XMIN) 10,10, 11 
Ll CALL OUTPUT 
mec 10 1000  ° °°. °&°» © eg ee ee 7 © 
memo o 12 1=1,NP 
meet) =X(1)/2.0— 
mr =0S5Q0(2)/A 


—-———- — -« - -_-—-— = -_-— = - oe - - ~ Se ee = - — 





mmm a EE eT RS ym a 


_——— se = ot ow 


a 


2 ee me ee ee -- = ‘ - _— oe 








_ WRITE Leg lOO Ne Cena ls ee ; . 
100 FORMAT (6E20.8) 
- WRITE (64101) CHIT ,NFORM,NMOVE,XTEST 
101 FORMAT(LPE2Z0.7/1 5X, TH NFORM= 1647X,7H NMOVE= 16, 7H XTEST= OPP. ee 
moo 10°13 a ev an — 
c MAKE PATTERN MOVES 
8 NMOVE=0 


18 NMOVE=NMOVE+L 

moe l4 [=1,NP 
eet, 1)=P(2,1)+DPC Ld ext 1) 
Mere) ( 1)=QSQUAR(1,1,1) 

mvs (2)=-0S9(1))13—e13—15 
Men e5Q(2)=050(1) 

woe 16 I=1,NP 
Mom (2, 1)=P(1,1) 

CALL ICLOCK(LPTIME) 
MeerrFiLPTIME.GT.ITIME) GO TO Yl 
17 GO TO 18 

END 


= —_— - in 
ee —s—“— — — a _ 


eS ee ee ee ee + ee - ee ee — —_— ——— ———— 





wee ~<a ~ cee ~~ eee ee ee =-— — ~ _ a _— — a 
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BOREL. «6 
matlS SUBROUTINE CALCULATES THE CORRELATION MATRIX FOR A LEAST = 
SQUARE ANALYSIS. nae 1) a — gi 


~ SUBROUTINE COREL 
DIMENSION A(5,5),;D0ER(5,20),AINV( 50,5)4B(50,1) 

~ DIMENSTON €(100),R(100),V0100) ,P(3,20) -X(20) -QSQ(3),DP(20),S(1500) © 
OIMENSTON GIVE(300) 





— COMMON GIVE,SrpEepReVe Pe Xs QSQ,NEsNPeIQ,JQ,KQ,ITOP,XTEST»sNFORM,NMOVE | 
COMMON /COMTIM/ITIME 
—NK=NEW _— 


DO 5 J=1,NP 

iy) > OK = 1g NK ee 
NE=K 

een ICLOCK(LPTIME) ~~ 
IF(LPTIME.GIFeITIME) GO TO 90 

SmeERCIT,K)=DERIVIT,K)Y °°» © 
mo10 l1=1,NP 

mee 10 J=1,f 

PINVil,J)=0.0 | 

— ete k= la NK || eee eee 
AINV( IT, J)=AINVO LT, J) +0ERV I,K) #®OER(IKIS(IR(K)#R(K) ) 

MeN TINUG™ 


lo eo oe we me re et ee ee ee et ee ee wee me ee ee ee ee ee ee ee ee ee ee ee i ee 


ue 


SemeeNP—) em SS 
—_— = Ly FM LE 
fee PP i=1+1 OU... OL 


DO ll J=IPLyNP 

evi lsJ)=BINVOJ, 1) 

Li CUNTINUE ~ a 

WRITE(6, 200) ((CAINV(1,J)_, L=leNP), J=1,NP) 

~ 200 FORMAT(36H INVERSE OF BEST CORRELATION MATRIX//(5E20.8)) 
on vee MATINV(AINV,NP,8620,0ETERM) 

meer i2 t=1,NP° © 

moel2 J=)1,NP 

meray |, J)=SQORTCAINV(I,J)) — 
Beebe CONTINUE 

- mime (6,L00)( (ACI sJ)y LTH=TRINP), J=l,NP) = ae 

100 FORMAT(44H BEST CORRELATION MATRIX (SQRT OF ELEMENTS)//(5€20.8)) 


me ee mr ge 


—S— ee ee ee ee ee 


rn emer rr rm rn rm mR tga pA cm 


90 RETURN ~ 
END 
oe ee ee 
me we ee cet weg ee ee ee ee ee nee a ee ee see ee me ee ee et ee ee eee ee ey ee ee ee ee ee ee ee on ee ee ee we ee ee —_ = We wee ee we ee ei ee ae ae 
ee me meee ee a ee ewe ee ee — ee me ee tte te ee em ee ee me we <r eee ee ee ee we oe oe ee ee eh at A Se agi ee ee a ee ee ee we ee ee - _- 
fete me PPP PA eee eal 
—_— = —— — ee ~~ -_ = — = — _——— — —_ a a am ee <2 = ow =0 <—_ _— a 
es -_ - = ——— _ ~ —~ _ — — — - 
a Ee a alee ae a —— ae ee ne ee 
_—- ~~ —_ —_ _— -_ — _ ee ~ —_— _ 


Tm ey ee ee ee Om oe re ee eg ee ee ee ee ee ee re ce ee ee —— = — —_— = — — es ee ee 





aA an on ae ee ae ee te es eee ee ee er ee eo eee ee eg ee me 


mm en me me en A A RN RM am ne ee ar a ee ee ———— 


Mery PUT... 
THIS SUBROUTINE CONTROLS THE OUTPUT PHASE OF SEARCH 
SUBROUTINE OUTPUT | 
~ DIMENSIUN E(100),R(100),V( 100) ,P (3,20) 5X(20) .QSQ(3),DP(20),8(1500) 
DIMENSION GIVE(300) 
— COMMON GIVEsSyEaRoeVeP 9X2 QSQeNE sNP_y IQ, 5SQ¢KQ,10P,XTEST sNFORM»NMOVE SA 
Sr =QS0(2)/A 





Meee XTEST=XTEST/2.0~ a. ; 7 
WRITE (6,200) 
OO FORMAT(38H FINAL VALUES OF PARAMETERS AND CHI ##2)00 ©0000 
WRITE (6,800)(P(2,1)+T=1,NP) h 


~ 800 FORMAT (6E20.8) 
WRITE (6,801) CHI »NFORMsNMOVE,XTEST 
~ 801 FORMAT(IPE20.7/15X_, 7H NFORM= I6¢7X,7H NMOVE= 16,7H XTEST= OPF9.1/) 
WRITE (6,300) (X(I), I=1l,NP) 
OO FORMAT(22H FINAL VALUES OF STEPS/(6E20.8)) °° & »° » ”»}» }©»™[W[V—— 
abl COREL 








ee ee ete ee — re - - ~ me a re ce ee rm re ee re re er ee ee es ee ee ee ee ee ee ee ee eee eee i a _ 
CALL THEORY(1 210 
9 9 9 
RETURN 
ee ee ee FE ee ee ere em ct ee ey ee ie ce oe ne ee ee ce ee et ee ee ee ee ee es es ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ~_~—— — — —_—— — 
NO 
me ee ee ee re re ee re me ce ce ee ec ee ere a ee en ee ee ee ee ew re re ee re ee ee we ee ee ee ee ee ee eee eee ee Oe -_-— a 
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THIS FUNCTION CALCULATES THE PARTIAL DERIVATIVE OF THE OUTPUT 


OF THE SUBROUTINE THEORY WITH RESPECT © TO SCHEER PARAMETER IP Al THE 
POINT JP. 


_DERIV3 [S FOR THE SPECIAL CASE OF THE SHIELDING PARAMETERS 
ONLY IN THE CASE OF C2 SYMMETRY. 
FUNCTION DERIV(IP,JP) 
~~ COMMON/COMDER/DOR(300),05(300) 
IF CIP-1)10,10,20 ; 
MERIT V=DR(JP) 
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C PUD T « << 
C JPLOT, THIS SUBROUTINE PLOTS A SINGLE VALUED FUNCTION 
© ON THE OFF-LINE PRINTER. THE ORDINATES ARE PLOTTED ACROSS THE PAGE _ 
C AND THE ABSCISSAS ARE PLUTTED DOWN THE PAGE 
= ARRAY TO BE PLOTTED (ONE DIMENSIONAL) 
NP= NO. OF ELEMENTS IN S 
NL= NO. OF LINES THAT PLOT IS TO COVER 
—  XO= LOWER LIMIT OF ORDINATES 
Xl= UPPER LIMIT OF ORDINATES 
[UNITS PER LINE e(SCAEE FACTOR). ~ 
DD= AN ALPHAMERIC TITLE FOR PLOT», 1LOAG 
LAB= NO. OF -THE PLOT, MUST START WITH 1 
WHEN LAB=1 THE PLOT SYMBOLS ARE READ. THE PLOT SYMBOLS, AA(1), 
~ CONSIST OF 7 WORDS THAT CONTAIN 6 BLANKS OR 5 BLANKS ANDO ONE | 
eieus SIGN. 
SeeSUBROUTINE JPLOT{S,NP,NL,XOyX19X—X1900;E AB) 
DIMENSION S(1500),GIVE(300),AA(9),AB(20),0D(10) | 
COMMON GIVE  § , § 
IF {LAB8-1)1000,1000,1001 
MEAD) «6(5,999) TAACI YS Tal eS) 0 eee 
999 FORMAT (9A6) 
wor 2D=NP/NL 
X2=(X1+X0)/2.0 
Me WRITE (6,980)0D,X0,X2,Xl © ‘<a. } [ 
980 FORMAT(18HL RESULTS OF JPLOT//1LOAG//L2X FT 2p 43K oF Te 24 43Xe Fle 2/ 
— Seo, 102H #@1ILTLLILIle LI File riiyiriiielriililireliliriliil l= 
POeePSITLILLITII# ILI lilielIIlilliqlelIIIifiiite) 
~ 00 998 J=1l,NP,JUD 
20 X*FLOAT( J) 4X1 _ 
Meee 94 [=15;17 } }§=©=© 2 2} 2 EE ESS Se 7 
994 AB{( I J=A AACT) 
IF (S(S)-X0)43944944 | 
43 AB(1)=AA(8) 
Mert 2)=AA(9) 
GO TG 990 
44 1F(S(S)-X1)46,46,45. 
45 AB(16)=AA(8) 
Beret t7)=AA(9) 
GO TO 990 
46 PS=100.0#(S(J)—-X0)/(X1-X0)+0.5 
PeS=PS 
i EST=0 
DO 997 [T=1,17 
“TFCIPS-1ITEST)1952,995,952 
eee FC IPS~6-I1TEST)950,991,991 __ 
~ 991 ITEST=ITEST+6 | 
997 CONTINUE 
S50 IDEL=I1PS-ITEST 
SOelO (1s2, 39400) LOEE 
995 AB(I)=AA(1) 
GO TO 990 
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met il }=AA(5) 

feo, TO 990 

ABLI)=AAL4) 

160 TO 990 

JAB TI )=AA3) 
mo 10 990 °° ~ 
ABT )=AAL2) 
iSO 10 990 
= AA6) 
meLTe (6,9897ABC,(ABIK),K=1l,1¢) °° 

FORMAT (7X)F6. alin Ee I 17A6) 
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§ CONTINUE 
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~ HIGHER J MAY BE HANDLED IF THE APPROPRIATE DIMENSION STATEMENTS | 
___ ARE CHANGED. THE ELEMENTS ARE CALCULATED FROM GIVEN CEF PARAMET. 


meee F M A T saaeas 
Pits SUBROUTINE CALCULATES THE CEF MATRIX ELEMENTS OF A RARE 
EARTH TON IN NTERM TERMS OF THE OPTICAL SPECTRA. YHE MAX J=15/2. 


INPUT DATA 
P(2,1)=CEF PARAMETERS IN THE FOLLOWING ORDER C20.C22,040,042, 


— 04-2 4C44404-47C60 9062 C6-24C64406-490661C6-6 


NTERM=NO. OF OPTICAL LSJ TERMS 
NJ=NO. OF DIFFERENT J VALUES IN TERMS a * 
XJ=MIN~e: J VALUE - 1.0 


Se ees a es ee ee. - - ~ . — - ~ 
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699° 
1 
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101 
999 








1000 
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102 
997 


See G96 F lige AX 
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C 
me )6hté<i~i~*S:SCKTTA:«WSJ:«CVALUE FOR ALPHA, BETA, GAMMA RESPECTIVELY, I AND J 
E 


S moetl,J.-K)=0. 


[W=DUMMY VARIABLE 
IDEG=1l2FOR INTEGRAL J, 2 FOR HALF INTEGRAL J 
THETA(N»M)=REDUCED MATRIX ELEMENTS, N TS ORDERYM IS TERM 
[IDENTIFICATION 
AJ(T)=J VALUE FOR THE ITH TERM  ——-__——————<“<—CSCSCSCSCS~—~— Sr 
NOP=NO. OF CARDS WITH OPERATOR EQUIVALENT MATRIX ELEMENTS FOR ~ 
_— EACH J VALUE DIVIDED BY 3 

OA, OBy OC(K,1,J)=1J OPERATOR EQUIVALENT MATRIX ELEMENT FOR 
GO FROM 1 TO 2J¢ly 1 IS INDEX FOR MJ=-J AND 2J+1 IS INDEX 
Oe es) = ae “- ae a 7. -<———o 
SUBROUTINE CEFMAT(1Q,JQ,KQ) 
DIMENSTON OA(8,16516),08(8,16116),0C(8416,16),THETA( 3,20) 
OIMENSION P(3,20), ANM (34442) 2AJ(20) »SIGN(2) 


~~ COMMON/COMTH/CEF(13,13,2) 


READ(5,699)(P( 2s), 1=1, 14) 
FORMAT(6E12.5) ~ Si — 
READ (5,100)NTERMsNJyXd¢ 1 Wy IDEG 
FORMAT (212,F10.55212) 

DU 999 M=1,NTERM 

READ (5,101) (THETA(N,M) »N=1,3),AN(M) 


FORMAT (3£20.5,F 10-5) 


+ te ee ——— = — a oe we 


seONTINUE rr 


DO 1000 1=1,8 

mer lo00 J=1,16 

wee1000 K=1,16 

fier l,.J,K)=0. 

il ,;Je,K)=0. 

CONTINUE 

DO 998 K=1,NJ 

READ (5,103)NOP 
FORMAT (12) 

D0 997 L=1l,NOP 

ae = (55102) Paar URN Ky 1 POOR CR pl) UCUIR TI) 
FORMAT (212,3F20.5) 
CONTINUE 

AK=K 
LMAX=AK4+XJ+0.52 
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uve = 
MIN=L¢l 
MeAX=2.08(AK+XJ)+1.01 
MAX=JIMAX41-—-L 
~~ DUO 996 N=MIN MAX ~ ia ii = ae 
NA=JMAX4+1—N 
DA(KsMAX,N)J=OA(KeNAoL) 7 
OB(K»MAX,N)=OB(KyNAgL) 
week ~MAXsN)=OC(K NA we 
996 CONTINUE 
Memee995 (=)],JMAX ©» © 0ST 
oe 995 N=l.L 
| Mee, N,L)=OA(K,LeN) © ee) ae . - 
OB(K»N,»,L)=OB(KeL,yN) 


mek ,N,L)=OC(K,LeN) 
995 CONTINUE 
meo08 CONTINUE =~” 
C SETTING UP OF ANM 
me DO 1001 I=1,3 
meelOO0l J=1.4 
meeeeoo) K=l,2  §° eee eee ge 0 = ag aed en 
Prom (1,J,K)=0. 


A a AY ET A S-Series ee ee ee 
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~~ ee re ee es ee ee eee ae ee we ee le -— ~~ ee —s- = Hie + = —- 


— eee ee oe 


_— 2 1A=0 Lk ot et eo _ i 
IS=1 
Mee 94 Nabe dS re ‘one 


~ MAX=N¢4)]_ 

DO 994 M=1,MAX 
e994 t=L,2 
Mett-N) 111,111,994 
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melt? CONTINUE oe a : —e 

| IF (M-L)994,993,993 

m993 [A=I1A+1 — eee ee cs a 
IF (LA-1Q)992,991,992 


992 ANM(N~M,L) =P(KQ, 1A) 
So 10 994 
mea ANM(N,4,L)=P(JQ,1A) 
994 CONTINUE 
— CALCULATE CEF MATRIX ELEMENTS _ 
| __SIGN(1)=1.0 
| perGN (2°)=-1.0 
if 1 
Me AJS(1)-xJ+.01 
—JMAX=2.08AJS( 1) 41.01 © 


mR ct et ee cee Rn RR RE ESS ARSE SEES 


me PO 903 J=l,JMAX ; as 

DO 903 K=1,JSMAX 
ee D0 903 t=1,2 ee) or Ong ee _- 
Per (o,Kh,ly=0.0 
903 CONTINUE : ——_. : 

DO 901 J=1,JMAX 
we 901l- Kaigd 

Beet JK) 80078017800 Eee 


me em te ree em ct cee ce te er ce mt en eR ter ory a A ER eee ce ee ee ee ne ee ee ee = re ee eee 
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——- al — a eee ——) 
801 M=l 
. GO TO 810 : 
800 IF (J-K-2)901,802,803 
—oent=2 © Oo _~ ON 
GO TO 810 
~~ 8C3 IF(J-K-4)901,804,805 ; ; 
804 M=3 
| GO TO 810 
805 I[F(J-K-6)901,806,901 
™ BCO6 M=4 ee | hae aa ee es a. ae 


810 CONTINUE 
- mo 3899 L=1,2 
CEF(JsKeLI=HSIGN(ILI*(THETACL, FL) *#O0ACT Ist eK) OANM(1LSM,LI*THETA(2,1)% 
- FL OB(TSsJ eK) #ANM(2sM,LItTHETAC3,T)*OC(TJsJrK)#ANM(3,MyL)) | 
CEF(K,JyLI=SIGN(L) #CEF(J,KyL) 
meg? CONTINUE a i 
901 CONTINUE 
-_ WRITE(6,1002) (P(2,;1), [=l,15) 
meer ORMAT(32HI VALUES OF CEF PARAMETERS USED/ BH A2OR2=F 10.3; TH eAge 
—  — IR2=F10.3,7H B22R2=F10.3, 7H AGORG=F10.3,7H AG2R4=FLO.3/8H B42R4= © 
foe LO 
he Se TH ASG4R4=F10.0347H B44R4=F10.23,7H AEORG=F10.3,7H AG2RE=F1O.37 
38H B&2RG6=F10.357H Ab4GRO=F10.397H BO4R6=F10.397H ASGRE=F10.3, 


_— oe oe ow ee ee 
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ce ee we ee ee oe -e 


47H B66R6=F10.3) 


RETURN 
pe me sen —_— E ~~ - = . ~~ sewe . —— - -_ me mee ee -_s 8 -= =— aw = ee ee ee _ _— - 
NO 
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fete OR OV Gs eo ecce ec oeneme 
THIS SUBROUTINE DIAGONALIZES CEF MATRIX PROVIDED BY CEFMAT AND 
CALCULATES(DELTA E)T, THE TEMPERATURE DEPENDENT QUAD SPLITTING 
~~ IN TERMS OF SHIELDING PARAMETERS. 
INPUT DATA 
NPOINT=NO. OF POINTS THAT ARE TO BE CALCULATED FOR PLOTTING 
IF NPOINT IS GREATER THAN OR EQUAL TO 150 THE TEMPERATURE 
DEPENDENT QUAD SPLITTING IS ALSO PLOTTED ON THE MOSELY 
PLOTTER 
SCALE=SCALE FACTOR FOR TEMPERATURE SCALE 7 UUW 
C20,C22=CEF PARAMETERS 
DO=LABEL FOR PLOTS 
YMIN»s YMAX=MIN. AND MAX. VALUES OF QUAD SCALE IN CM/SEC 
~~ SUBROUTINE THEORY(1Q, JORKG, OP) 
COMMON GEV&SS,EeReVeP eX: QSQ,NEr NP TQeJQeKQ, TOP yXTEST »NFORMyNMOVE GA 
~COMMON/COMDER/DR(300),0S(300) ~— ars sar rr 
COMMON/COMHER/EVENUP »W(13) 2H(78),SPARE(79) 
“COMMON /COMTH/CEF (13,1352) , ; ; 
DIMENSION Q22(13),QX¥(13),D0(10) »C(13413)EX(13,100) 
— DIMENSION GIVE(300) © lel 
DIMENSION E(100)5R(100),V( 100) »P (3920) X(20) »OSQ(3)50P (20) 1911500) 
DIMENSION EVECT(13),EVECR(13)  °&°+&2&=&7} ee 
COMPLEX H,C,CTEMP,CMPLX 
MemmPLleEX CSQRT, CONJG, CABS 
IF (NFORM)1,2,1 
-_ Meer CEFMAT(1,2,2) 
oe 3 [=1,13 
MeL =CEFCI,I~1) ° ©=— 00000 
3 CONTINUE 
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DO 4 K=JP1,13 
Mer =CMPLX(CEF(U,GK,1L),CEFL IR Ke2))  °@ @&«&«°«#«8©§©)—ClC)”6)C)h)©r?.”D)6WVWVU™U US 
+] 
4 CONTINUE 
SPARE(79)=0.0 
DFLTA=1.0E-10 
moet HERM{(H, 13,C,0,D0ELTA,I1T) 
“CC . .©6MULTIPLY C BY PROPER PHASE FACTOR 
mee 1000 1=1,13 
MEMPZREAL(C(1,7)) 
[F(TEMP.NE.0O.0) GO TO LOOL 
Mr eEMP=CSORT(-CONJG(C(I,12))/C(I,2)) 
fm 10 1002 
meeoimeCTEMP=CONJG(C(I,7))/CABS(C(1,7)) ° 
Moe2 DU 1003 J=1,13 
1003 C(l,J)=CTEMP#C(1,J) 
1000 CONTINUE 
C REARRANGE ENERGIES AND WAVE FUNCTIONS 
HOS] = 42 
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Pe il=[+} 
DO 5 = Pues 
[F(WCl)-W(5))5,5,7 
Meet EMP=W(I)~" °° &# # # - Sn ae ee 
W( IT) =EW(J) 
~ WIJ) =TEMP 
mos K=1,13 
MerTeEMP=C{t1,K) 
mets K)=C(5,K) 
fc{(J,K)=CTEMP 
5 CONTINUE 
—_— RENORMALIZE THE ENERGIES | 
TEMP=W(1) 
mo 9 [=1,13 © 
9 W(1)=wWwllI)-TEMP 
C- ~=3DU|W CALCULATE MATRIX ELEMENTS OF 35Z*##2-J(J+1) ANDO 372(JPeae2¢+JMen2) 
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ee er ee ee oe 


mee iO J=1,13 
Q22(1)=Q2Z(T) #REAL(CONIG(C (Td) ) *C (Ty S)) #03208 (FLOAT (S)=“7.0) 
l #®*£2-42.0) 


a we 68 om ae Oe 88 oe ae 





10 CONTINUE ——- 
fee) f=1,13 
MY (1) =O .0 0 EE 
mo 11 J=3,13 
— °° AJ=J-7 . OS a aes ee eee 
rey 1 )=QOXV( 1) +1. S*REAL (CONIG(C(1,J)}*Cli,J-2)+C(1,J)8 | 
~ 1} CONJG(C(I,J—-2))) *SQRT((6.0+AJ)#(5.04AJ) #(7.0-AJ) #(8.0-AJ) ) 
, 11 CONTINUE = - ——————t—te 
mmr0(S,100) NPOINTsSCALEsC20,C22 26242. °° 
mega FORMAT (14,F10.5/2F20.5) 
—- B=0.695056 SS eee Se ne” re 
E20=0. 769854 
—  HCE=0.8610308E-5 _ 
R2=0.19 
Mee eR 3=75.5 De eee 
ALPHA=1.0196651E-02 
Me SCALCULATE THE PARTITION FUNCTION °°” °° ©&©&« © ©™©)0DCUCCUSU SS 
1 IF (10P)12,13,12 
Mrs MI=NPOINT 2 = TS 
, Meet 
ei? Ml=NE I i | 
IF (NFORM)15,14,15 
> 14 DO 19 M=1,Ml _—. e a | dA a - : ‘ Sa : 
16 T=V(M) 
Sees DO 19 N=1,13 I 
EX(NyM)=EXP(—-W(N)/(B#T)) 
meeeeir9o CONTINUE | Pa ae <a ae . 
DIRL=0.5#E2Q#4.08#HCE#C20/R2"(-1.0) 
; DIRZ2=0.5#E20#4.08#HCE#C22/R2"(-1.0) 
__C4F=0. S#E2QSALPHA#R3#(-1.0) oe 
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15 CONTINUE 
FEED IN THE PARAMETERS AND CALCULATE 
GAMMAR=P(KQ,1) 
me 10 (50,51),1Q 

~~ §0 GAMMAR=P({JQ,_,1) 

) GO TO 52 

mo, GAMMAS=P (JQ, 2) 
52 CONTINUE 

men 20 J=1,M1l — 

| PetlOP)70,%],70 

— T1 T=SCALE*FEGAT(Y) 
DO ee l=l, 3 

72 EX(1,LI=AEXP(-W(1)/(B8#T)). 
wep =} 

meer, IF (LOP.NE.0) JAP=Y 

aR 1=0.0 
oa 
21 PART=PART+EX(1,JAP 


ee eet ee ee ee oe 





a ih oe al oo a oe ~-—— = 


ee 
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ee Ne ee Ppp ne ae nee eee re ee 
ee 


Sx YA=0.0 
Peto) 22 I=1,13 
M7 n=QZ2ZA+0ZZ(1)*#EX( 1, JAP) 
22 QXYA=QXYA+QXY(T)#EX(1,3AP). 
SA=C4F #Q2ZA/PART 
~ SAA=SA*#(1.0-GAMMAR)+DIR1L*GAMMAS 
IF (DIR2)2001,2000,2001 
Meeeooes(j)=SAA ©}. 
mo TO 2002 
B 2001 SB=C4F2QXYA/PART | 
SBB=SB#(1.0-GAMMAR)+DIR2*GAMMAS 
S(J)=SORT(SAA##241.0/3-08SBBe#2) © 
moe 1F(10P)23,20,23 
23 CONTINUE | 
met J )=-(SAA*SA+SBB#SB/3.0)/5 
DS (J)=(SAA#DIR1+SBB#DIR2/3.0)/5 
20 CONTINUE 
IF(10P)24,25,24 
25 WRITE(6,200) IT 
~~ 200 FORMAT ( 32+ 
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i 48HO ENERGY 

e11?7H (CM-1) MJ= aaa =o 
3 O a +2 +3 

me) 30 l=1,13 - 


DON 31 J=1,13 
EVECR( J) =REAL(C(1,J)) 
EVECI( J) =AIMAG(C(19J)) 
WRITE(6,201L)W(1),(EVECR( J) s 
FORMAT(F10.3,6X 13F8.3/16X 13F8.3) 
CONTINUE | 
WRITE(6,202)QZ22,QKY 


31 


201 
30 
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NUMBER OF I[TERATIONS IN HERM 


J=1,13),(EVECI(K), 


em er rr 


THE QUAD SPLITTING 
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Se 
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wl oe 
MIXING COEFFICIENTS? 
7 3 -2 -1 
+4 _ +5 +6// \ 
K=1,13) 
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2 FORMAT {(18HO 
Bet J-*4#2)//(6E20.8) } 
~READ(5,203) ODyYMINs YMAX — all 
Bee ORMAT (LOAG/2F20.5) gm ee 
WRITE(6,204) DD,SCALE,(S(1),1=1,M1) 
04 FORMAT(22H1 RESULTS OF PLOT OF 10A6/10X 15H SCALE FACTOR= 
mee lO.2//(2X 20F6.3)) 
BerrF(NPUINT.LT-150) GO TO 32 
CALL CPLOT(SeM1yM1,0D,YMAX,YMINe 1) 
2 CALL JPLOT(S,M1l~-120,YMIN,YMAX, SCALE, 0.0,00,1) 
24 RETURN 
No 
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HITEMP 


The required card decks are: 


THE OR Ye 

DERIV a2 

SEARCH 

QSQUAR see listings for QTAVE 
OUT Pus 

COREL 


JP LoOe 


The listings of decks 1-2 follow. 
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C THO ELOURTY 5 

Cc THIS SUBROUTINE CALCULATES THE QUADRUPOLE SPLITTING IN THE FORM 

C OF A/T+h 

oe |6™~*~*é<~«wWNPUT COAT AC” - es | ee i 
Ee TEMP=INITIAL TEMPERATURE WHERE PLOT IS TO START 

———_ SCALE=MAX. TEMPERATURE/40 . 7 

c DD=LABEL FOR PLOT 
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Se 
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SUBROUTINE THEORY(1Q,JQ,KQ,I10P) 


COMMON GIVEr Se Er RaVaPy Xe QSQsNEaNPe IQs JQeKQe1OPeXTEST yNFORM, NMOVEs A 
—COMMON/COMOER/S4FeSLAT = 


DIMENSION GIVE(300),S4F(0120) 
DIMENSION 00D(10) 


DIMENSTON E(100),R0100)5V0100) sP (3420) eX(20) sQSQ(3),0P (20) 4S5(300) 


~~ ee —_— 


IF {(NFOGRM)1_1¢2 

Il R2=0.19 

~ R3=74.0° 00° 
ALPHA=1.0201LE-2 


_ — — _ eee ee ee ne ee ee ee ee ee ne ee ewe ee ea es Kee - _ 


ALPHAZ=ALPHA®ALPHA | 
J=6 


feo —Je(J+lLie(2eJ+lLoe(2eJalye(2eStoy/ 5 Lee 


RM1L=1.0 
B=0.695056 ~ 
E2Q=0. 769854 
mer-0.8610308E~-5 


Oe ad ee ee re _ — 


ee ee eee oa ~ 


ANUM=-0. DFE2QeRZ#RZeALPHAZ*SUM#RM1/( 13. 0*B) 


-_— ee me ee eee ee ee ne ee 


ALAT=0.5#E2Q*#HCE#4.0 
READ(5,100)TEMP,SCALE 


100 FORMAT (2F20.5) 


oO oe 


oe ee ee me 


2 ANM=P(KQ,1) 
GO TO (344),1Q 

MD j0,1)0 
GO TO 5 


Meee GAMMA=P(JQ,2) 00° 7 0 


— +s tee 


—-— ~~ — - 


5 SLAT=~-ALAT 
IF (1LOP)6,7,8 

¢ M1=4) 
so 10 9 

8 ML=NE 

moo 109 

6 Ml=l1 

moO 10 M=1,M1 
Mev lLOP)11,12,13 


8 ee eee ee es _— —- - ee 


Mie t=V (NE) CoC eaal | ogee 


GO TO 14 


Meee 1 =TEMP+SCALE*FLOAT(M) = = = = ——— 


mee CON LUE 


GQ TO 14 

13 T=V(M) 

14 S4F(M) =-ANUM/T 
S(M)=-ANM#S4F(M)-GAMMA#SLAT 


MEM A=P(KO,2)°> 000 ee . » 2 |. 


i eR a Qa | cer 





mr 10OP)15,16,15 
meme 2, 2)=P(2,2)/P(2,1) 
WRITE( 6,200) P25) ) 9 Gece) 
90 FORMAT(28HO SQRT(A204"#2 + 1/3A22082)= F10O.2,22H (1-GAMMA)/(1- 
1GMA)= FlO,cw 
READ(5y101)00 ——— 
LOL FORMAT (10A6) 
WRITE(6,102)00, TEMP, SCALE ae 
02 FORMAT(50HL PLOT OF TEMPERATURE VS QUADRUPOLE SPLITTING IN LOA6, 
“140HO NO. OF POINTS @=741°>~~ INITIAL TEMP" = F10c1, 7 
Bae SCALE FACTOR & F10-1) aa 
Berna (6,103) (S(LI@ ee L Mie eS a ee 
3 FORMAT(2X 20F6.3) a 
~~ CALL JPLOT(S441441,0.0,10-0,SCALE,TEMPyDD, 1) © 7 7 
15 RETURN 
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MERIV 2.-. 
THIS FUNCTION CALCULATES THE PARTIAL DERIVATIVE OF THE JUTPUT 
OF THE SUBROUTINE THEORY WITH RESPECT TO THE PARAMETER IP AT TH! 
POINT uP. 
OERIV2 IS FOR THE SPECIAL CASE OF THE SHIELDING PARAMETERS 
ONLY. 
FUNCTION DERIV( IP, JP) 
DIMENSION E(100),R(100)_ V6 100) »P (3920) 2X20) »QSQ(3) 51390) 
DIMENSION GIVE(300) 
DIMENSION S4F(120) 
COMMON GIVEs,SpEeReV_9 Pe Xe QSQ_NE NP, IQ, JQ,KQ,10P,XTEST »>NFORM,NMOV: 
COMMON/COMOER/S4FySLAT 
Ment THEORY(1,22251) 
IF (1P-1)10,10,20 
10 DERIV=-S4F (JP) 
ico 10"30— 
20 DERIV=-SLAT 
30 RETURN 
__ END 
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